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in  content;  2)  the  biopsy  specimens  were  re-read  using  a  new  set  of  standards  so 
■that  corrections  for  thermally  induced  tissue  shrinkage  could  be  calculated;  3) 
graphical  representations  oi  the  data  contained  in  the  revised  data  base  were  made 
as  a  help  in  understanding  the  relationships  embodied  in  the  data;'  4)  an  empirical 
multidiscriminant  model  was  written  which  predicts  either  gross  (clinical)  grade 
or  histopathologic  grade  given  parameters  such  as  heat  flux,  exposure  time,  skin 
temperature ,  and  the  like;  5)  an  analytical  model  was  developed  to  circumvent 
the  problems  associated  with  the  multidiscriminant  model  regarding  the  inability 
to  use  a  flux-time  profile;  6)  optimization  of  the  analytical  model  was  accomplishec 
using  data  in  the  data  base,  intraskin  time-temperature  profiles  and  such  data  as 
could  be  found  in  the  literature,  e.g.  University  of  Rochester,  Moritz  and 
Henriques  and  Stoll,  as  performance  criteria.  Abstracts  of  10  additional  reports 
covering  various  aspects  of  the  project  in  more  detail  and  a  listing  of  BRNSIM, 
the  analytical  model  are  appended..'  Conclusions  and  Recommendations  include: 

1)  The  USAARL  Porcine  Bum  Data  Base  is  ready  for  release. 

2)  Data  handling  and  graphing  programs  are  available  to  facilitate  further  study. 

3)  BRNSIM,  which  predicts  bum  depths  over  a  wide  range  with  reasonable 
accuracy,  should  be  tested  against  the  Aerotherm  manikin  model. 

4)  Skin  cooling  by  blood  and  tissue  water  boiling  should  be  further  studied. 

5)  BRNSIM  should  be  tested  using  Stoll's  extensive  human  bum  data  base. 
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1.0 


INTRODUCTION 


This  is  the  final  report  of  a  data  analysis  and  model  development  study 
based  on  information  collected  in  USAARL's  Thermal  Analysis  Program.  (1,  2) 

The  overall  technical  objective  was  to  provide  managers  with  a  clinically 
valid  algorithm  for  evaluating  thermally  protective  fabrics  with  poUt-nilal  for 
military  use.  The  basic  question  which  managers  need  answered  i  -.  whether  a 
particular  fabric  will,  when  exposed  to  an  accurate  simulation  of  a  nost-crash 
fire,  reduce  thermal  transfer  to  clinically  safe  levels,  and,  if  r--  . '  •  how  long? 

Colateral  to  this  basic  question  is  whether  the  fabric  itself  add;-  hazard  by  becoming 
a  secondary  thermal  source.  If  a  particular  fabric  becomes  a.  secondary  source,  it 
is  then  important  to  know  the  degree  to  which  it  increases  the  hazard. 

In  order  to  answer  these  questions,  a  logical  approach  would  be  to  expose 
fabrics  to  a  controlled  thermal  source  which  accurately  simulates  a  post-crash  fire 
*nd  to  measure  the  thermal  transfer  through  the  fabrics  with  skin  simulants  or 
heat  flux  transducers.  This  approach,  if  carefully  carried  out,  would  provide 
accurate,  reproducible  engineering  data  on  thermal  transfer  through  the  fabrics. 
However,  it  would  not,  by  itself,  provide  clinically  meaningful  data.  Calories 
per  second  transferred  is  not  biologically  significant  information.  What  is  signifi¬ 
cant  is  the  resulting  tissue  damage.  It  is  this  thermally  induced  damage  and  not 
the  heat  transfer  per  se  which  incapacitates  or  kills  the  army  aviator . 


1.1  EXPERIMENTAL  RATIONALE 

The  original  Thermal  Analysis  Program  (1,  2)  was  set  up  to  determine  the 
relationships  which  exist  between  thermal  transfer  and  injury  production.  The 
experimental  rationale  was  as  follows;  Heat  flux  transducers  and  the  skin  of 
pigs,  as  a  human  skin  analog,  were  exposed  to  the  sams  accurately  simulated 
post-crash  fire  produced  by  a  specially  designed  shuttered  furnace  burning  JP-4 
fuel  (3).  Both  sensors  and  skin  were  exposed  either  bare,  or  while  covered  oy 
layers  of  test  fabric.  Thus  the  sensors  gave  an  accurate  indication  of  heat  flux  as 
a  function  of  time  with  direct  exposure  or  as  modified  by  fabric  layers.  The  skin 
tissue  samples  gave  a  valid  measurement  of  clinical  damage. 


1.2  DATA  COLLECTION 

July.  1972  to  January,  1973  a  large  body  of  experimental  data  was 
collected  using  the  USAARL  porcine  bioassay  technique  (1,  2,  3)  relating  thermal 
flux  and  exposure  time  to  bum  depth.  Physical  thermal  sensors  were  also  ex¬ 
posed  to  the  simulated  post-crash  fires  in  hopes  of  establishing  a  correlation  be¬ 
tween  their  output  and  the  bursts  received  by  the  pigs. 


1 . 3  ANALYSIS  BY  UTRI 


During  the  collection  period  and  shortly  thereafter,  ITT  Research  Institute, 
under  contract,  analyzed  the  data  that  was  being  collected  and  developed  a  mathe¬ 
matical  model  for  predicting  bum  injury  (4)  based  on  earlier  work  by  Henriques 
(5) ,  Stoll  and  Greene  (6) ,  Weaver  and  Stoll  (7) ,  Mehta  and  Wong  (8)  and  Morse 
et  al  (9).  UTRI's  model  could  accept  either  a  constant  heat  flux  as  an  input 
parameter  or  a  table  of  heat  fluxes  as  would  be  generated  from  the  heat  flux 
sensor  either  in  a  laboratory  test  or  as  part  of  an  instrumented  manikin  used  in 
"fire  pit"  testing.  But,  during  the  data  collection  and  model  building,  both  ItTRI 
and  USAARL  personnel  observed  that  the  data  was  noisy  and  that  the  model  could 
be  improved.  There  were  also  some  items  of  data  either  not  available  during  this 
period  or  not  utilized  by  UTRI  such  as:  1)  intraskin  temperature  recordings;  2) 
heat  flux  values  derived  from  analog  output  of  the  heat  flux  sensors  which  had 
been  recorded  on  FM  magnetic  tape;  3)  bum  depth  readings  which  allowed  for 
corrections  to  be  made  accounting  for  the  shrinkage  observed  to  occur  in  the  more 
severely  burned  tissues;  4)  additional  measurements  of  tissue  water  and  density. 

1.4  WORKING  HYPOTHESIS 

The  hypothesis  of  this  study  was:  given  the  output  of  physical  thermal 
sensors  and  the  actual  bums  experienced  by  pigs,  it  should  be  possible  to  develop 
a  mathematical  model  which  can  convert  the  output  of  the  sensors  to  predicted 
bum  depths  in  reasonable  agreement  with  the  observed  bum  depths.  To  the 
extent  that  this  were  possible,  it  would  be  reasonable  tc  use  such  a  physical 
thermal  sensor  and  the  mathematical  model  thus  derived  as  part  of  a  screening 
methodology  for  thermally  protective  fabrics. 

1.5  OBJECTIVES 

The  specific  objectives  accomplished  toward  developing  a  clinically  valid  fabric 
screening  method  are  as  follows:  1)  the  previously  collected  data  base  was  reviewed 
for  error  both  in  coding  and  in  content;  2)  the  biopsy  specimens  were  re-read 
using  a  new  set  of  standards  so  that  corrections  for  thermally  induced  dssue 
shrinkage  could  be  calculated;  3)  graphical  representations  of  the  data  contained 
in  the  revised  data  base  were  made  as  a  help  in  understanding  the  relationships 
embodied  in  the  data;  4)  an  empirical,  multi-discriminant,  model  was  written 
which  predicts  either  gross  (clinical)  grade  or  histopathologic  grade  given  parameters 
such  as  heat  flux,  exposure  time,  skin  temperature,  and  the  like.  However,  since 
this  model  required  either  a  constant  input  flux  or  the  use  of  total  flux,  its  use 
is  limited,  making  the  model  incapable  of  using  the  thermal  exposures  experienced  at 
the  "fire  pit"  test  facility  of  Natick  labs;  5)  an  analytical  model  was  developed  to 
circumvent  the  problems  associated  with  the  muitidiscriminate  model  regarding  the 
inability  to  use  a  flux- time  profile;  6)  optimization  of  the  analytical  model  was 
accomplished  using  data  in  the  data  base,  intraskin  time-temperature  profiles  and 
such  data  as  could  be  found  in  the  literature,  e.g.  University  of  Rochester, 

Moritz  and  Henriques  and  Stoll,  as  performance  criteria. 
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1.6  SUMMARY 

In  summary  then,  the  objectives  were  to  validate  the  data  base  collected  by 
the  staff  at  USAARL,  to  look  at  the  data  in  graphic  presentation,  to  ascertain 
functional  relationships  and  problems  with  the  data  and  to  develop,  ultimately, 
a  model  which  would  take  the  output  from  heat  flux  sensors  and  accurately  predict 
the  bum  depth  which  would  have  occurred  had  the  sensor  itself  been  replaced  by 
living  skin ,  The  remainder  of  the  report  will-  deal  with  the  work  accomplished  to 
meet  each  of  these  objectives  taken  in  the  order  just  presented. 


2.0  DATA  BASE  VALIDATION 


Data  from  the  thermal  project  were  recorded  in  three  forms;  1)  handwritten 
records,  e.g.  gross  (clinical)  grading  sheets;  2)  stripchart  recordings  of  temperature, 
heat  flux  and  shutter  system  state  as  a  function  of  time;  3)  analog  voltages  repre¬ 
senting  temperature  and  heat  flux  recorded  on  FM  magnetic  tape.  At  the  start  of 
this  contract  the  written  records  had  been  coded  on  computer  oriented  record  forms 
and  key-punched  on  computer  cards.  These  cards  were  read  into  data  files  which 
then  became  the  data  base.  None  of  the  data  on  FM  magnetic  tape  had  been  utilized 
although  it  had  been  digitized  and  stored  on  digital  magnetic  tape  as  part  of  the 
contract  performed  by  UTRI  (4) . 

In  order  to  form  a  data  base  containing  pertinent  information  from  all  sources , 
which  was  reasonably  free  of  error,  the  following  tasks  were  accomplished. 

2.1  ELIMINATE  KEYPUNCH/TRANSCRIPTION  ERRORS 

The  entire  data  base  was  read  item  by  item  and  compared  with  the  original 
data  records  (clinical  grading  forms,  pathology  reports,  etc.).  Errors  were  eliminated 
and  those  observations  which  did  not  match,  e.g.  a  gross  grade  of  15  with  a  burn 
depth  of  2C0|im,  were  further  checked  against  the  photographic  records  and  biopsy 
material  to  resolve  the  mismatch. 


2.2  RE-READ  BIOPSY  MATERIAL  TO  ALLOW  FOR  SHRINKAGE  CORRECTION 

It  became  apparent  late  in  the  data  collection  period  (10)  that  biopsies 
from  deep  dermal  burns  showed  a  marked  thermally  induced  shrinkage  while 
more  superficial  burns  showed  some  edema  and  swelling.  The  edema  and  swell¬ 
ing  was  slight,  approximately  3-5%,  and  not  as  dramatic  as  seen  in  human  burns 
(11).  The  shrinkage,  on  the  other  hand,  was  at  least  40%  and  went  as  high 
as  216%.  Clearly  correction  for  such  a  large  distortion  of  measured  burn  depth 
would  enhance  the  chances  of  a  model's  ability  to  predict  bum  depth. 

With  this  in  mind,  the  entire  set  of  more  than  1600  biopsy  slides  was  re-reaa 
using  descriptive  criteria  as  before  (1)  but  adding  several  new  depth  measurements 
as  shown  in  Figure  1.  By  recording  both  normal  skin  thickness,  (A  +  B),  and 
burned  skin  thickness,  (E) ,  as  well  as  "depth  of  bum"  C  it  was  possible  to  cal¬ 
culate  a  "corrected  bum  depth"  as  follows:  (A  +  B)  -  C  [(A+B)/El  «  C.B.D. 

Measurements  A,  B,  C,  D,  E  and  the  corrected  bum  depth  were  added  to 
the  data  base1 .  The  histopathological  grading  definitions  along  with  mean  experi¬ 
mentally  observed  bum  depths  are  found  in  Table  I.  A  second  transformation  was 
made  to  the  bum  depths  as  follows:  [Corrected  BD/(A+B)]  *  2000  «  N.B.D.  This 

1 A  more  detailed  discussion  of  the  histopathology  grading  system  can  be  found  in 
USAARL  Report  78-11,  June,  78.  Abstract  App.  B. 
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THRESHOLD  BURN  DEPTH 


A  NORMAL  EPIDERMAL  THICKNESS 
8  NORMAL  DERMAL  THICKNESS 
C  •BURN  DEPTH' 

0  DERMAL  THICKNESS  AT  BURN 
E  TOTAL  SKIN  THICKNESS  AT  BURN 


Figure  1.  Schematic  Diagram  of  Porcine  Skin  Showing  Location 
Measurements  A  through  E  recorded  for  each  biopsy  specimen. 


"normalized  burn  depth"  corrected  for  the  large  variability  in  total  skin  thickness 
(A+B).  Normalized  burn  depth  was  added  to  the  data  base. 

2.3  CORRECTIONS  TO  ORIGINAL  DATA 
2.3.1  Shutter  Dynamics. 

There  were  four  different  shutter  propulsion  systems  employed  during  this 
study.  In  a  study  of  the  fire  simulator/ shutter  system  (see  Abstract,  Appendix  B) 
correction  factors  were  derived  for  each  shutter  system  and  each  bum  site  location. 
These  correction  factors  were  applied  to  each  exposure  time  in  the  data  base  90 
that  the  time  represents  the  best  estimate  of  actual  exposure  time  ignoring  compli¬ 
cating  factors  such  as  day  to  day  variations  in  exposure  time  caused  by  changes  in 
the  rolling/ sliding  friction  of  the  shutter. 


2.3.2  Initial  Skin  Temperature. 

Initial  skin  temperature  has  been  shown  to  be  an  important  factor  in  determining 
the  energy  requirements  for  a  given  level  of  burn  damage  (12) .  Skin  temperatures 
were  recorded  at  one  of  the  bum  sites  for  each  exposure  beginning  5  September 
1972.  No  such  recording  was  made  prior  to  that  date.  A  multiple  linear  regression 
analysis  (see  Table  II)  was  used  to  derive  a  relationship  between  environmental 
factors  (order  of  bum,  time  of  day,  temperature  at  Carins  AAF,  %  humidity,  wind 
velocity  and  %  cloud  cover)  and  measured  skin  temperature  for  pigs  burned  sub¬ 
sequent  to  5  Sept  1972 .  This  relationship  was  then  used  to  calculate  a  predicted 
skin  temperature  for  those  pigs  burned  during  July  and  Augu3t.  This  temperature 
was  inserted  into  the  data  base. 


2.3.3  Hair  Removal. 

Two  different  hair  removal  procedures  were  used  during  the  course  of  the 
experiment  and  this  resulted  in  mean  hair  length  of  3.2mm  for  pigs  burned  in  July 
and  August  and  1.75mm  for  the  remaining  pigs.  In  the  first  case,  only  about  93.75% 
of  the  incident  radiation  gets  to  the  skin  while  in  the  latter  case,  96.5%  does  so  (4). 
Hair  length  tor  each  pig  is  in  the  data  base  to  aid  in  correcting  the  incident  heat 
flux. 


2.4  DATA  MANAGEMENT  COMPUTER  PROGRAMS 

In  order  to  edit,  retrieve,  sort,  list  and  otherwise  manage  the  data  in  the 
data  base,  and  to  retrieve  and  format  data  from  digitized  tapes,  e.g.  records  of 
sensor  output  digitized  from  FM  magnetic  tape  recordings,  some  12  programs  were 
written,  primarily  in  FORTRAN  IV.  Their  complete  description  including  source 
listings  and  sample  output  may  be  found  in  a  report  entitled  "Thermal  Analysis 
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Program:  Computer  Programs  for  Data  Base  Manipulation"  (see  Appendix  5  for  an 
abstract) .  A  one  page  sample  listing  of  the  primary  data  base  is  shown  as  Tab  Li 
QX.  The  group  (GBP)  number  in  the  last  column  is  a  quality  control  rating  from 
1  to  9.  A  rating  of  4  indicates  prime  quality  data  which  was  collected  during  the 
svunmer  months  when  the  experimental  protocol  was  undergoing  considerable  develop¬ 
ment.  For  a  complete  listing  of  the  data  base  which  includes  a  discussion  of  all 
the  quality  scores,  see  a  report  entitled  "Experimental  porcine  bum  injury  data  base: 
listings  and  graphic  representation."  (see  abstract  in  Appendix  B). 

3.0  GRAPHIC  PRESENTATION  OF  THE  DATA  BASE 

A  computer  graphing  program,  PLOTS,  allows  the  us<ar  to  plot  any  two  variables 
on  arithmetic  or  log  axes  parameterized  by  a  third  variable  if  desired.  The  follow¬ 
ing  examples  give  an  indication  of  complexities  and  the  basic  variability  inherent 
in  the  data.  Verticle  bars  indicate  sample  means  and  the  individual  data  points 
are  indicated  with  a  symbol  instead  of  a  bar  for  ±1  standard  deviation  or  ±1 
standard  error  of  the  mean.  The  first  four  plots  show  gross  (clinical)  and  micro 
(histopathologic)  grades  versus  normalized  and  corrected  bum  depth.  In  all  cases 
there  is  a  logical  progression,  though  not  linear,  from  low  grades  and  shallow  depth 
to  high  grades  and  deep  bum  depths. 

Both  systems  apparently  lack  precision  in  distinguishing  degrees  of  shallow, 
i.e.  epidermal  bum.  It  turns  out  that  this  is  mostly  a  function  of  the  bum  depth 
measurement  scheme  in  which  C  (see  Figure  1)  is  measured  from  fat/dermal  border  up 
to  epidermal  border  and  stops  there.  Thus  there  is  no  really  accurate,  quantitative 
depth  measurement  when  damage  is  limited  to  some  fraction  of  the  epidermis.  For  more 
discussion  of  this  and  other  "Pitfalls  in  the  use  of  the  USAARL  Histopathology  Grading 
System",  see  abstract  in  Appendix  B. 

The  next  two  plots  relate  bum  grade  to  total  heat  flux.  Again  there  is  a 
logical  progression  with  more  severe  bums  requiring  more  energy  to  produce  them. 
Neither  of  these  representations  reveals  the  difference  between  the  bums  resulting 
from  one  second  at  3  cal/cm1*sec“1  versus  8  seconds  at  1  cal/ cm1 /sec.  Since 
reciprocity  has  been  shown  not  to  hold  (13),  there  is  reason  to  believe  that  the 
picture  is  more  complex.  Moreover,  none  of  these  plots  takes  into  account  such 
uncontrolled  variables  as  initial  skin  temperature,  length  of  time  anesthetized, 
length  of  hair,  etc.  All  these  other  variables  have,  most  probably,  contributed  to 
the  noisiness  of  the  data. 

The  next  plot  compares  the  gross  and  micro  grading  systems.  The  correlation 
between  them  is  quite  good  (r1  3  .99)  from  gross  grade  6  through  16  and  micro 
grade  4.S  to  10.  There  are  so  few  data  points  below  this  level  that  it  is  not 
possible  to  draw  any  conclusion  regarding  the  relationship  between  gross  and  micro 
grades. 

The  next  four  plots  are  an  attempt  to  explore  the  relationship  among  three 
variables ,  flux,  exposure  rime,  and  bum  grade  or  depth.  The  log  exposure  time 
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vs.  log  flux  parameterized  by  bum  grade  plots  were  suggested  by  Alice  Stoll  (14) 
based  on  her  experience  with  pain  thresholds  and  threshold  blister  bums .  In 
her  plots  these  Iwo  parameters  form  straight  lines  which  are  essentially  parallel. 
The  USAARL  data  suggests  a  shift  to  the  upper  right  as  a  function  of  burn  depth 
as  might  be  expected  but  there  also  appears  to  be  a  slight  change  in  slope  of 
the  "best"  (drawn  by  eye)  straight  line  from  one  burn  level  to  the  next. 

The  last  plot  (figure  12)  is  a  different  view  of  the  same  data,  namely:  nor¬ 
malized  depth  versus  total  flux  parameterized  by  exposure  time.  While  there  is  a 
progression  to  the  data,  it  is  largely  obscured  by  the  internal  variability. 
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4.0  EMPIRICAL  MODEL 


The  lack  of  really  clear  cut  relationships  in  the  foregoing  graphical  analysis, 
prompted  the  development  of  a  multidiscriminant  model  which  takes  into  account  not 
only  time,  flux,  total  flux  and  burn  grade,  but  includes-  furnace  wall  temperature, 
skin  temperature  and  skin  color  as  well.  The  results  of  applying  this  model  to  pre¬ 
dict  gross  burn  grade  is  shown  in  Table  IV.  The  generalized  Mahalanobis  D-square, 
the  equivalent  of  a  Chi  square  for  this  sort  of  analysis,  is  highly  significant  showing 
that  it  is  possible  to  predict  with  reasonable  certainty  the  gross  bum  grade  knowing 
exposure  time,  flux,  furnace  wall  temperature,  total  flux,  initial  skin  temperature, 
and  skin  color. 


Multidiscriminant  analysis  was  also  applied  to  predict  micro  grades  -with  a  some¬ 
what  lower  D-square  but  generally  similar  results.  Both  models  were  used  to  flag 
outlying  data  points  for  further  checking  against  original  records.  However,  neither 
model  lends  itself  to  use  with  an  instrumented  manikin  and  fire  pit  testing  where  the 
time-heat  flux  profiles  recorded  by  the  data  acquisition  system  are  not  constant  but 
change  as  the  overlying  fabrics  react  to  the  fire  which  is  itself  variable In  this  case 
an  analytical  model  is  more  suitable. 

5.0  ANALYTICAL  MODEL 

Several  years  ago  Weaver  and  Stoll  (7)  proposed  an  extension  of  Stoll's 
earlier  model  (6)  to  heat  fluxes  higher  than  used  in  obtaining  the  experimental  data 
upon  which  the  earlier  model  had  been  based.  They  also  found  that  the  effective 
conductivity  changed  during  the  exposure  and  subsequent  cooldown  period.  Takata 
(15),  using  preliminary  data  from  USAARL's  Thermal  Project  (the  uncorrected  version 
of  the  current  data  base) ,  formulated  a  model  which  not  only  predicted  threshold  burns 
but  deep  bums  and  tissue  water  boiling  as  well.  Building  on  the  work  of  Henriques 
(5) ,  Stoll  and  Greene  (6) ,  Weaver  and  Stoll  (7) ,  Mehta  and  Wong  (8)  and  Takata 
(15)  we  formulated  an  analytical  model  as  follows: 

For  thermal  exposures  of  interest,  skin  is  essentially  opaque  to  thermal  radia¬ 
tion  and  can  be  considered  to  transfer  energy  internally  by  conduction  only,  since 
exposure  durations  are  no  longer  than  the  minimum  response  times  reported  for 
increased  thermoregulatory  system  activity  (16).  Consequently,  thermal  energy 
transfer  in  skin  can  be  described  by  the  heat  conduction  or  Fourier  equation.  In 
rectangular  coordinates,  the  Fourier  equation  may  be  written  as  follows: 


P 

r_  9T  _  9  fv  9 Tv  , 

Cp  aT  -  Si  «  ^  1 

where, 

P 

= 

density,  gm/cm5 

Cp 

= 

heat  capacity,  cal/gm/°C 

K 

S 

thermal  conductivity,  cal/cm-sec.°C 

(1) 


T  =  temperature,  °C 

x  =  distance,  cm 

q  =  energy  source,  for  the  first  nodal  volume,  cal/cm*/sec 

Since  skin  is  considered  to  be  opaque  to  radiant  energy,  and  since  the 
source  term  is  due  only  to  radiant  energy1 ,  equation  (1)  applies  only  to  the 
surface  of  the  skin.  For  all  conditions  in  which  x  >  0,  equation  (1)  reduces  to 
the  following: 

§T  _  8  .  dT 

p  P  8t  ”  8x  K  8x  ^ 

Solution  of  equations  (1)  and  (2)  requires  two  boundary  conditions  for  x, 
preferably  at  x  —  0  and  x  —  L ,  and  initial  conditions  at  t  —  0  for  all  positions 
0<X<L.  If  one  assumes  that  there  is  no  backward  flux  of  thermal  energy  at 
x  =  0  (all  conduction  is  into  the  skin),  then  the  energy  flux  at ’x  s  0  is  zero  and, 
consequently,  8T/8X  =  0.  Similarly,  if  the  problem  assumes  that  an  adiabatic 
back  well  condition,  prevails  at  X  =  L,  the  fatty  tissue,  then  the  net  flux  out  of 
the  system  at  X  -  L  is_  .0,  or  9T/9X  =  0.  These  two  boundary  conditions  indicate 
that  the  system  is  closed  and  that  all  thermal  energy  added  to  the  system ,  (XX<L , 
is  distributed  within  the  system  and  cannot  escape.  J  ‘ '  Js‘  w  ” 

Initial  conditions  are  established  by  specifying  a  uniform  temperature  for  all 
locations,  (KX<L,  at  time,  t  =  0. 


Consequently,  the  system  may  be  defined  by  the  following  mathematical  model: 


rn  21  3  i  rv  9Ts  . 

P  CP  3t  3*  1 


@  x  =  0 


p  p  aT  8x  CKa7 


@  0<X>L 


T  =  To,  (KX>L,  t  =  0 


Initial  conditions 


~  =  0.  x  =  0,  0<t>x 

8x  - 


Boundary  conditions  1 


- —  =  0,  x  =  L,  0<t>x 
8x  - 


Boundary  conditions  2 


^  A  simplifying  assumption  based  on  the  predominance  of  the  radiate  mode  of  heating 
May  be  less  valid  with  fabrics.  In  actuality  a  correction  is  made  to  q  to  account 
for  convective  heating,  surface  absorptivity,  and  attenuation  of  radiant  heating  by 


6.0  SOLUTION  OF  MATHEMATICAL  MODEL 


An  analytical  solution  to  equation  set  (3)  was  not  considered  feasible  due 
to  the  variable  nature  of  q,  Cp  and  K;  so  explicit  differencing  methods  of  numerical 
analysis  were  employed  to  solve  the  equations.  Several  investigators  working  with 
linear  systems  have  found  that  the  Crank-Nicholson  six  point  implicit  differencing 
method  provided  an  excellent  numerical  solution  (17).  For  the  solution  of  equation 
set  (3),  the  mathematical  model,  it  was  decided  to  apply  the  Crank-Nicholson  method 
to  the  second  order  partial  derivatives  and  corresponding  explicit  methods  to  the 
Erst  order  partials. 

The  grid  work  in  Figure  13  is  a  representation  of  the  differenced  system  from 
X  =  0  to  X  =  L  (J's)  and  t  =  0  to  t  =  y  (i's). 

The  Crank-Nicholson  technique  involves  averaging  the  value  of  the  dependent 
variable  bvir  the" i  and  i  *  1  row  at  a  constant  j  position.  The  second  order  deriv¬ 
ative  is  thfen  evaluated rat ‘the  (j ,  i  +  1/2)  position.  A  forward  difference  formulation 
is  applied  tb  the'  3T/3t  term;  to  match  the  same  position.  .7  '.7 

~  »  :a  Q 1. 1 .  5  «  *  t!  :  C 

Trhe  'ab'ov^  'described  implicit  differencing  method'  ii'nof&fft$ *ihe  characteris¬ 
tics  or  stability- 'khd  convergence.  Correct  increment  sii’eS  yfeld- reliable  convergence. 
The  mcdel^ah-  implemented  in  FORTRAN  IV  vising  solutitSii  ^Chhic^ies  of  Thomas  as 
described  by  Bruce  et  al  (18).  ,y"e  2 

This  initial  model  was  revised  to  allow  energy  flux  across  the  surface,  x  =  0, 
during  heating,  convective  heat  loss  at  the  skin  surface  during  cooling  and  heat 
transfer  into  deep  tissues  including  conduction  into  fat,  convective  cooling  via  the 
blood,  tissue  water  boiling,  a  temperature  gradient  from  surface  to  fat  and  a  gradient 
of  thermal  properties  based  on  measured  tissue  water.  For  a  complete  description 
of  the  model  including  source  listings,  sample  output  and  users  manual,  see  a 
USAARL  report  entitled  "BRNSIM  -  An  Analytical  Model  for  Predicting  Thermal 
Cutaneous  Injury",  abstracted  in  Appendix  B.  The  model,  BRNSIM,  is  run  inter¬ 
actively  with  the  following  variables  changeable  for  each  run: 


TABLE  V 


MODEL  PARAMETERS  CHANGEABLE  INTERACTIVELY 


Initial  Surface  Temperature  (TEMPIC3) 


Density  (p) 

Thermal  Conductivity  £cr  each  Node  (BK(i)) 

i  • 

» 

Node  Depth  (BL  =  0.22  cm) 

Heat  Capajdty  for  each  Node  (Cp^) 
IteCiEcn  lntirval  (AK  -  0.01  sec) 

f 

j 

Number  of  Nodes  (JINC  =  12) 


Exposure  Time  (ETIME) 


Total  Time  (ITIME) 


Water  Boiling  Temperature  (°C) 

Blood  -  Control  Factor  for  Convective  Tissue  Cooling  by  Blood 
TEMP  B  -  Difference  Between  Back-wall  Tissue  Temperature  and  TEMP1Q 
Absorptivity  of  Surface 
Incident  Flux,  Q,  (cal/ cm* /sec) 

Nxtra  Nodes  -  Number  of  Interpolated  Nodes  Between  Surface  and  Node  2 
Used  for  superficial  bums. 


Damage  Rate  Constants  PL1,  PLN1,  PL2,  PLN2,  DEI,  DE2 
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From  the  relationship  for  first  order  kinetics  assumed  to  apply  in  damaging 
tissue  protein  we  have:  (5) 

total  (^TIME 

damage  =  J  ^/dt  «.  \  dfl/dt 
/  0  i  ETIME 


=  PL  +  PLN  -  DE*  1 

(T+273) 

Thus  for  damage  calculations  the  following  constants  are  entered:  (15) 

PLX  (44°C  -  50°C)  =»  1.46  PL2  (50°C  -  100°C)  ’»  2.24 

PLNX  (44?C  -  50°C)  *  147.37  PLN2  (50°C  -  100°C)  =»  239.4-7  "r  -'•  - 

DEX  (44°C  -  3Q°C;)  =  50,000  DE2  (50°C  -  1O0»C)  =  80,000 

on?  4 

The  program  outputs  dft/dt,  for  each  node  at  each  time  step,  total  is. 
damage  for  each  node  and  a  threshold  depth,  where  £1  =  1.  This  depth,  found 
using  inverse  interpolation  on  two  or  three  £l's  nearest  1  using  either  y  or  log(y). 

Since  its  first  presentations  (19,  20)  BRN’SIM  has  undergone  further1  develop¬ 
ment. 


damage  rate  =  =  pe“AE/®-T. 


dt 


if  P  =  N  x  10?  and  AE/R  =  DE 


then  In  tt*  =  InN  +  yin  10  -  •  i 

R  T 


dt 


7.0  THERMAL  PROPERTIES  OF  SKIN 

Measurements  of  the  water  content  of  pig  skin  as  a  function  of  thickness  were 
made  on  split  thickness  skin  samples  from  several  pigs. 

Given  a  table  of  measured  values  of  water  content  as  a  function  of  skin 
thickness,  a  least-squares  cubic  polynomial  was  fit  to  the  data  and  water  content 
as  a  function  of  depth  was  computed  from  the  formula: 

w(T  -  d)  -  I  (\VT  -  WT-d)  ♦  WT-d 


where  T  is  the  total  thickness  of  a  slab,  W-j.  is  the  fraction  of  water  computed 
from  the  cubic  equation,  d  is  the  thickness  of  a  thin  slab  at  a  depth  T-d,  and 
Wt-d  is  the  traction  of  water  above  the  thin  slab. 

Thermal  properties  of  the  tissue  were  computed  from  the  equations:  (21) 
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1)  density: 


Y  a 


^w 

yw 


W. 


f  wp 

Yf  YP 


-1 


2)  heat  capacity:  Cp  =  Ww  Cp^  +  Wf  Cpf  +  Wp  Cp? 

3)  thermal  conductivity:  K  -  y  ^w  +  kf  Wf  +  kp  Wp 

Y™  yf  YP 


where  the  subscripts  w,  f,  and  p  refer  to  water,  fat,  and  protein,  respectively. 

Wn  is  the  mass  fraction,  yn  the  density,  CPn  the  heat  capacity,  and  kn  the  thermal 
conductivity  of  the  respective  components.  Values  of  the  various  terms  used  were: 


yw  =  1  gm/cc  *'•  "cPw'"  1  cal/SPn~°C 
yf  -  0.815  gm/cc  .»  ^  . Cp^  =  0.55  cal/ gm-°C 
yp  =  1.54  gm/cc.  c -  Cg^*.=  0.26  cal/gm-°C 


k^sIil^SXlO-3  cal/cm-sec-°C 
kf  *  4,5X10"4  cal/ cm-sec-°C 

Or  .  a  r  0  'jc  - 

j£pr=  .4t32y.gr4  cal/cm-sec-0C 


Fat  and  protein  were  assumed  to  be  present  is 


a-:  ..  us.?.  ■;  - 


,  r..:C*K.- 
Y  •  i  ■ 


i  (i  -  Ww), 


'"'C1'  0  s  s-  r  '"V5;'  ■'  .. 

equal  amounts  so  that": 

'b'SLb  axy-cruo  cangvir. 
s  biu  ebon  riot? 
no  noiaaiofHs*'?su  e?v 


and  the  resultant  equations  were: 


y  =  (6.18277X1G-*  Ww  +  .938 172) ~x 

K  =  y(l. 08432X10“®  Ww  +  4.15684X10*4) 
Cp  =  .595  Ww  +  .405 


Using  the  equations  above,  the  profile  of  thermal  properties  was  calculated 
for  skin  depths  of  from  81  to  2290 jiin.  A  linear  extrapolation  of  tissue  water 
content  from  a  depth  of  81}im  to  the  skin  surface  was  made  using  a  stratum  comeum 
watar  content  calculated  from  Rushmer  et  al  (22)  and  the  ambient  %  humidity  during 
the  experimental  phase  of  the  project.  This  c.dculated  water  profile  was  used  to 
complete  the  calculation  of  thermal  properties  profile  from  81}im  to  the  skin  surface. 
These  new  thermal  properties  replaced  those  chosen  by  Morse  et  ai  (9)  and 
used  during  previously  reported  simulations  (19,  20).  See  the  report  entitled 
"Thermal  Properties  Calculated  from  Measured  Water  Content  as  a  Function  of 
Depth  in  Porcine  Skin"  as  abstracted  in  Appendix  B  for  additional  details. 


8.0  INTRASKIN  TEMPERATURES 

In  earlier  simulations  (19,  20)  it  became  apparent  that  unless  the  temperature 
calculations  reasonably  represented  what  actually  occurred  in  the  skin,  adjustment 
of  the  values  for  PL,  PLN  and  DE  in  the  damage  equation  to  match  a  few  data 
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points  would  not  be  likely  to  res-—  in  a  model  which  works  well  for  all  cases. 
Fortunately  eleven  intraskin  temperature  profiles  were  recorded  on  FM  magnetic 
tape.  These  voltage  records  were  digitized  and  converted  to  tables  of  temperatures 
at  100  samples  per  second.  Figure  14  present?  the  one  page  summary  report  from 
a  simulation  of  the  exposure  of  Pig  294RF  to  a  3.47  cal  cm"1  sec-1  fire  for  3.02 
seconds.  Note  that  boiling  occurred  (confirmed  by  blister  formation.  Figure  15)  and 
that  the  surface  reached  a  maximum  of  128.724°C.  Predicted  threshold  depth  was 
1520|im.  Three  observed  temperature  profiles  are  overlayed  on  the  calculated 
temperature  profiles  (for  nodal  depths  of  0,  220,  440. . . .22QQ^m)  in  Figures  16,  17 
18.  The  oscillations  in  the  observed  temperature  profile  are  most  probably  due 
to  a  "hunting"  in  the  autoregulation  of  tissue  perfusion  by  blood.  The  frequency, 
for  example,  is  similar  to  that  seen  in  studies  of  microcirculation  (23). 

The  next  series  of  figures,  19-22,  shows  a  simulation  in  which  the  intraskin 
temperature  gradient,  cooling  by  blood  and  water  boiling  are  turned  off.  This  time 
the  temperature  profile  does  not  fit.  The  threshold  depth  is  similar  only  because  the 
heat  flux  is  higher.  (The  comparable  simulation  with  gradient,  blood  and  boiling 
turned  on  had  a  Max.  Ten^».  »  130.945,  Threshold  depth  *  1585  and  a  ftnaT  time 
of  100).  The  next  simulation  (BRNSIM  3  vs  M12VB012.0,  see.. Figures,, 23- ?6)  is 
identical  with  the  previous  except  that  gradient,  blood  and  boiling  are^hnmed  on 
and  the  recalculation  of  skin  thermal  properties  subsequent  to  water  bqilipg^is 
inhibited.  Note  that  this  simulation  does  not  fit  well  either.  Thus/^ohly  fixe  complete 
model  reasonably  simulate^  fixe  temperature  profile.  y.  V  •  1 

:  « •  .*v  ■■ 

Table  VI  presents  a  comparison  on  burns  observed  in  the  bioassay  test  with 
predictions  of  Model  3  (BRNSIM) .  This  table  is  an  updated  version  of  one  pre¬ 
viously  published  (20)  and  includes  normalized  bum  depth  in  addition  to  corrected 
burn  depth.  The  effect  of  normalized  burn  depth  in  reducing  the  RMS  Error  can 
be  seen  in  Table  VII.  One  problem  with  Model  3  was  noted  in  running  case  296RF. 
The  normal  version  of  the  model,  labelled  M12VB0120  or  Model  3,  includes  a  recal¬ 
culation  of  thermal  properties  subsequent  to  boiling.  This  is  done  percepifiously 
following  the  cessation  of  boiling  and  in  cases  of  high  heat  flux  and  moderate  or 
longer  exposure  times;  it  acts  as  a  pulse  of  heat  flux  causing  an  instability  as  can 
be  seen  in  the  temperature  plot  in  Figures  27  and  28.  A  run  with  BRNSIM  3,  which 
is  Modal  3  with  the  recalculation  routine  inhibited,  results  in  a  prompt,  but 
perhaps  too  rapid,  cooling  of  the  skin  following  shutter  closure  (see  Figures  29  and 
30).  This  is  not  a  problem  in  those  cases  when  boiling  does  not  occur.  It  can 
be  fixed  by  redesigning  the  recalculation  algorithm. 


9.0  PREDICTIONS  FROM  SENSOR  DATA 

The  next  set  of  simulations  used  the  digitized  and  converted  output  from 
Aarotherm  Thermomaxx  Sensors,  a  file  of  flux  values  as  a  function  of  time,  to  drive 
the  model.  Table  Vm  summarizes  twenty  simulations  and  compares  the  predicted 
depths  with  both  normalized  and  corrected  depths  and  Figure  31  shows  the  results 
of  one  simulation.  Data  produced  by  the  older  version,  Model  II,  from  Knox  et 
al  (20),  is  also  included.  As  noted,  columns  and  Illg  were  run  at  slightly 


MODEL  HOME  OB  DESCRIPTIONS  MODEL  3  PIS  294RF  DB  VALUES  BBS*. 613  *  INCIDENT  FLUX  6  APB  79 


SKIN  DIFFUSION  DATA 
INPUT  PARAMETER  LIST 


TEMP  IS*  34.9768 
DENS*  1.00060 
01*  3.47660 

BL*  0.220060 
AK*  0.I06600E-01 
U  INC*  12 
TEMPB*  3.3600 
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BOIL*  100.130 
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MAXIMUM  TEMPERATURE  •  128.724 


THRESHOLD  DEPTH  •  1320. 

FINAL  TIME  «  180.00 


Figure  14 


Figure  15.  Intraskin  thermocouple  (0.003",  "located  superficially")  shown  prior 
to  burn  (left)  and  subsequent  to  exposure  to  3.47  cal'cm'^sec-1  for  3.02  seconds 
(right) . 

Gross  grade  =  13  New  Micro  grade  =  8  Threshold  depth  =  1465jjim 
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MODEL  NAME  OR  DESCRIPTION!  M12VS0128  PIG  294  *27.23.29  WITH  TEMPR.  BLOOD  *0  AND  BOIL-250  4  APR  79 


SKIN  DIFFUSION  DATA 
INPUT  PARAMETER  LIST 


TEMPIO* 

34.9780 

DENS* 

t. 80000 

01-  3 

.74000 

BL-  0.220000 

AK«  0. 

1000A0E-01 

JINC- 

12 

TEMPB* 

8.nnA0 

ABSORB- 

0.013000 

BOIL- 

230.000 

PL2- 

2.24 

PLM2* 

239.47 

PI.  1  - 

1.46 

PUII* 

147.37 

DE2*  08000.00 

DEI-  50000.00 

ETIME- 

3.02 

1TIME* 

108.00 

NXTRA- 

0 

BLCOO* 

0.0000 

FLUX  FILE  I.O.t  8.80  2 


FLUXCD* 

1  3.740  2  3.748 


H- 

0.891O4E+8I 

W- 

0.43214E+00 

W-* 

0.2S090E-0I 

p- 

0.72442E+61 

p- 

0.73770E+O1 

0* 

0.74955E+01 

14  • 

0 . 9S073E+22 

AT  DEPTH  UN  MICRONS)- 

0.1123351 

w  • 

0.2OI34E+17 

AT  DEPTH  UN  MICRONS)- 

200.000 

II  • 

0.4I470E+I2 

AT  DEPTH  UN  MICFCNS)  • 

400.000 

14  • 

0. 1 I413E+09 

AT  DEPTH  UN  MICRONS)  ■ 

600 . 000 

II  • 

O.32S22E+0S 

AT  DEPTH  UN  MICRONS)- 

000 . 000 

N  • 

0.30033E+04 

AT  DEPTH  UN  MICRONS)- 

1000.00 

14  « 

0. 18037E+03 

AT  DEPTH  UN  MICRONS)- 

1200.00 

N  « 

0.89IO4E+01 

AT  DEPTH  UN  MICRONS)- 

1400.00 

14  • 

O.46214E-00 

AT  DEPTH  UN  MICRONS)- 

1S00.00 

W  » 

O.2S890E-OI 

AT  DEPTH  ( IN  MICRONS)- 

1000.00 

14  • 

0.31424E-02 

AT  DEPTH  UN  MICRONS)  ■ 

2080.00 

14  • 

0. 00000 E +00 

AT  DEPTH  UN  MICRONS)- 

2288.00 

I1AXIKIM  TEMPERATURE  •  158.243 

THRESHOLD  DEPTH  •  1547. 

FINAL  TIW  •  38.12 


Figure  19. 


46 


•SECO 


MOTEL  HOME  0*  DESCRIPTION!  BRHSIFO  FIG  294  *27,29.29 


FLUX  FILE  9.09  2 


FLUXt I)* 

l  3.740  2  3.740 


u- 
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0.74955E+0I 
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UM 
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H  • 
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MAX1NJM  TEMPERATURE  •  ! 
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THRESHOLD  DEPTH  -  1453. 

FINAL  TINE  •  33.93 


OB  CORRECTED  VALUES  ABS-.613  VS  M12VMI29  5  OPR 

SKIII  DIFFUSIOM  DOTO 
INPUT  PARAMETER  LIST 


TEFne* 

34.9709 

DENS* 

1.80009 

01*  3.74000 

BL*  0. 

220000 

AK*  0. 

1000006*01 

J  INC* 

12 

TEMPB* 

3.3600 

ABSORB* 

8.613009 

80 IL* 

100. 138 

PL2* 

2.24 

7LH2- 

239.47 

PL1* 

1.46 

FUI1* 

147.37 

DE2*  80000.00 

DEI*  50009.00 

rriME* 

3.02 

ITIME* 

100.00 

NXTRA* 

0 

BLOOD* 

0.0010 

0. 1 12333E-06 
200.009 
400.000 
600.000 
80O.0OO 
1000.00 
1200.00 
1400.00 
1600.00 
1800.00 
2000.00 
2200.30 


Figure  23. 
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this  simulation  the  routine  to  recompute  thermal  properties  subsequent  to  water  boiling  was  turned  off  (see 
for  discussion)  . 
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Heat  Flux,  Q,  =  Q  incident  to  sensory  with  skin  absorptivity 


different  heat  fluxes.  The  most  nearly  correct  flux  would  be  incident  flux  with 
an  absorptivity  of  .613.  Once  again  the  RMS  error  generally  favors  normalized 
depths  and  would  be  even  better  with  proper  flux  input  (see  Table  IX) . 

Once  again  the  error  in  shallow  burn  depths  is  pointed  up .  The  depths 
for  SIM31213  should  be  compared  with  the  micro  grades  of  4  and  5  indicating 
from  partial  to  complete  epidermal-dermal  separation.  Using  Table  I,  this  would 
give  a  normalized  burn  depth  of  from  125  to  200pm,  much  more  in  line  with  a 
Model  mB  prediction  of  141pm. 


10.0  COMPARISON  WITH  WEAVER  AND  STOLL  (7) 

The  results  of  just  two  simulations  are  presented  here  in  Figures  32  and  33. 
First,  Model  3  (M12VB0120)  was  run  for  Stoll's  human  burn  data  (.4  cal  cm*1  sec-1. 
Absorptivity  =  0.94,  ETIME  =  5.6  sec  and  initial  skin  temperature  of  32.5°C) 
using  a  slightly  lower  than  normal  value  for  blood  (.0007  vs.  .001)  and  Stoll's 
integration  constants  (PL1,  PLN1,  etc.).  The  resulting  calculations  for  £1  at  8 Op 
are  different  by  9.4%  and  the  total  time  by  less  than  0.4%.  In  Model  3  the  blood 
cooling  is  turned  on  linearly  from  0  to  maximum  during  the  Erst  20  seconds. 

There  is  some  indication,  see  Figures  20,  21,  and  22,  that  the  cooling  by  blood  does 
not  visibly  take  effect  until  about  16  seconds  after  the  beginning  of  exposure. 

An  indication  that  a  slightly  different  model,  in  which  the  blood  is  turned  on 
from  37  to  45  seconds  after  the  beginning  of  exposure,  may  be  better,  is  shown  in 
Figures  34  and  35.  Here  the  predictions  are  different  from  Stoll's  by  0.8%  for  total 
damage,  £1,  at  80pm  and  by  0.35%  for  total  time.  This  is  very  good  agreement 
considering  the  differences  between  the  models  —  Stoll:  measured  human  skin 
surface  temperature,  no  blood  cooling,  thermal  properties  constant  as  a  function 
of  depth,  changes  in  conductivity  as  a  function  of  time  and  human  skin  thermal 
properties  versus  Model  3:  calculated  pig  skin  surface  temperature,  cooling  by 
blood,  pig  skin  thermal  properties  variable  as  a  function  of  depth,  and  initial 
thermal  gradient. 

Substitution  of  Takata's  (15)  values  for  the  activation  energy  (DE)  and  fre¬ 
quency  factor  (PL  and  PLN)  in  the  first  case  above  results  in  a  predicted  depth  of 
57.1pm  compared  with  105.2pm.  This  indicates  that  different  values  for  DE,  PL  and 
PLN  should  be  used  for  the  epidermal  and  dermal  nodes.  This  combination  of  factors 
has  not  been  tried. 


11.0  UNANALYZED  SENSOR  DATA 

One  set  of  experimental  data  collected  during  the  1972/73  period  remains  to 
be  analysed.  Both  Aerotherm  Sensors  and  Fabric  Research  Labs  skin  simulants 
were  subjected  to  the  heat  from  a  4  inch  thick  steel  plate  mounted  in  the  furnace. 
Unfortunately  the  only  record  of  these  experiments  is  in  the  form  of  strip  chart 
recordings.  The  peak  temperatures  of  the  sensors  were  analysed  by  ETRI  (4) 
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to  ascertain  their  accuracy  and  reproducibility.  When  these  recordings  are  man¬ 
ually  read  and  digitized  then  the  voltage-time  profiles  can  be  converted  to  flux-time 
protiles  and  used  to  drive  the  model.  It  will  be  important  to  complete  this  analysis 
since  some  of  the  sensors  were  covered  with  fabric  spaced  i"  away  from  the  surface, 
and  represent  the  only  such  data  available. 

12.0  DISCUSSION 

The  data  base  has  been  reviewed  and  corrected  and  is  now  ready  for 
release  to  those  who  wish  to  pursue  model  development  on  their  own.  As  can 
be  seen  in  the  graphic  presentations  of  the  data,  there  is  still  a  great  deal  of 
scatter  which,  to  some  extent,  is  unavoidable  in  experiments  of  this  type.  There 
are  statistical  tricks  for  dealing  with  scatter  but  the  emphasis  throughout  this 
project  has  been  to  attempt  to  understand  the  sources  of  the  observed  variability. 
Uncontrolled  variables  such  as  skin  temperature,  time  between  induction  of 
anesthesia  and  exposure,  disagreement  among  pathologists  regarding  gradings, 
and  unrecorded  variation  in  shutter  performance  all  have  had  an  effect  which  cannot 
be  removed. 

An  empirical,  multidiscriminant  model  was  developed  to  predict  clinical  or 
micro  grade;  but  further  development  of  it  or  some  more  appropriate  model  was 
temporarily  shelved  in  favor  of  analytical  model  development  when  it  was  realized 
that  the  immediate  need  was  for  a  model  compatible  with  both  laboratory  testing 
and  fire  pit  testing  of  thermally  protective  fabrics.  As  we  better  understand 
the  relationship  among  total  flux,  exposure  time  &nd  burn  depth,  a  simplified, 
empirically  based  model  may  become  available.  For  example,  if,  as  Stoll  suggested 
(14) ,  plots  of  log  exposure  time  versus  log  flux  parameterized  by  burn  grade/ 
depth  results  in  a  family  of  well  behaved  curves  or  straight  lines,  then  such  plots 
migit  prove  useful  in  screening  fabrics.  The  first  try  at  this  approach  was  not 
entirely  satisfactory.  The  data  base  and  plotting  routines  are  now  available  to 
investigate  this  more  fully. 

An  analytical  model  was  developed  which  takes  into  account  tissue  water 
boiling,  a  temperature  gradient,  a  thermal  properties  gradient,  skin  absorptivity , 
either  constant  or  variable  (tabulated)  heat  flux  corrected  for  attenuation  by  hair 
and  a  10%  convective  component.  It  has  been  shown  to  predict,  with  reasonable 
accuracy,  burns  in  pigs  when  they  are  exposed  bare  or  protected  by  experimental 
fabrics  and  when  the  heat  flux  is  a  constant  or  when  it  is  variable  and  derived 
from  output  of  an  Aerotherm  sensor, 

13.0  CONCLUSIONS  AND  RECOMMENDATIONS 

1.  The  USAARL  Porcine  Burn  Data  Base  is  ready  for  release. 

2.  Data  handling  and  graphing  programs  are  now  available  which  facilitate 
the  further  use  and  study  of  these  data. 


3.  A  multi-discriminant  model  exists  which  predicts  clinical  and  micro  burn 
grades  but  requires  further  development  and  study  to  be  useful. 

4.  An  analytical  model,  BRNSIM,  has  been  developed  which  predicts  burn 
depths  over  a  wide  range  with  reasonable  accuracy. 

5.  It  is  recommended  that  the  recently  detected  flaw  in  BRNSIM  be  corrected 
and  that  BRNSIM  be  further  tested,  not  only  with  data  from  the  USAARL  data  base 
and  the  University  of  Rochester,  but  also  with  data  from  the  instrumented  manikins 
as  used  in  die  Natick  Fire  pit. 

6.  It  is  recommended  that  BRNSIM  be  tested  against  the  Aerotherm  manikin 
model  in  its  ability  to  predict  burns  over  a  wide  range  of  conditions. 

7.  It  is  recommended  that  skin  cooling  by  blood  and  tissue  water  boiling 
be  further  studied  to  clarify  the  dynamics  of  these  phenomena. 

8.  It  is  recommended  that  BRNSIM  be  further  tested  using  Stoll's  extensive 
human  burn  data  base. 
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BIBLIOGRAPHY  -  THERMAL  ANALYSIS  PROGRAM 
Abstracts  of  Reports  and  Papers  in  Preparation: 

Knox,  F.  S.,  Peter  Saurmilch ,  S.  C.  Knapp,  Thomas  Wachtei,  George  McCahan  and 
Lynn  Alford.  "A  Fire  Simulator /Shutter  System  for  Testing  Protective  Fabrics  and 
Calibrating  Thermal  Sensors"  (USAARL  Report  79-4)  Ft.  Rucker,  Alabama:  U.  S. 
Army  Aeromedical  Research  Laboratory,  March,  1979. 

The  design,  construction,  calibration  and  use  of  a  JP-4  fuel,  shuttered  furnace 
is  described.  Based  on  a  NASA  design  this  furnace  simulates  the  radiative  and  con¬ 
vective  thermal  environment  of  a  postcrash  fire  in  rotary  winged  aircraft.  Heat 
fluxes  ranged  from  0.5  to  3.6  ±  3%  calories  per  square  centimeter  per  second  with 
steady- state  furnace  wall  temperatures  from  519°C  (967°F)  to  1353°C  (2450°F)  and  a 
radiative /total  flux  ratio  of  approximately  0.9.  A  pneumatically  propelled,  wai._L 
cooled  shutter,  mounted  in  a  rolling  animal  carrier,  controlled  the  exposure  of 
pigs  and  thermal  sensors  to  die  fire.  An  electronic  data  acquisition  and  control  sys¬ 
tem  is  also  described.  This  system  automatically  controlled  the  opening  and  closing 
of  the  shutter  and  provided  strip  chart  and  FM  magnetic  tape  records  of  exposure 
time,  furnace  wall  temperature,  heat  flux,  and  sensor  output.  Sources  of  error  in¬ 
cluding  nonuniformity  of  flame  front  and  shutter  dynamics  are  discussed.  Methods 
of  animal  handling,  burn  grading  and  photographic  documentation  are  introduced  along 
with  a  brio*  description  of  some  nine  experimental  protocols  carried  out  using  this 
fire  simulator  shutter  system. 


Knox,  F.  S.  HI,  Thomas  L.  Wachtei,  S.  C.  Knapp.  "Experimental  porcine  burn 
injury  data  base:  listings  and  graphic  representation." 

Pigs,  as  human  skin  analogs,  were  exposed  to  simulated  postcrash  fires  of 
various  intensitites  for  various  durations.  Some  pigs  were  exposed  bare.,  others 
with  blackened  skin  and  still  others  protected  by  thermally  protective  fabrics.  For 
each  exposure  the  location  of  the  burn,  exposure  time,  heat  flux,  initial  skin 
temperature,  date,  clinical  and  histopathological  evaluations  of  burn  injury,  cor¬ 
rected  and  normalized  bum  depth  and  a  quality  control  number  are  listed.  The 
data  base  contains  over  45,000  items  and  summarizes  most  of  the  porcine  burn  data 
collected  using  the  USAARL  fire  simulation  furnace/ shutter  system.  The  report  also 
includes  graphic  representations  of  items  within  the  data  base  such  as  log  exposure 
time  versus  log  flux  parameterized  by  gross  grade,  gross  grade  versus  normalized 
burn  depth,  new  micro  grade  versus  normalized  burn  depth,  gross  grade  versus 
corrected  burn  depth,  new  micro  grade  versus  corrected  burn  depth,  gross  grade 
versus  new  micro  grade,  new  micro  grade  versus  total  flux,  gross  grade  versus 
total  flux,  normalized  burn  depth  versus  total  flux  parameterized  by  exposure  time. 
These  plots  will  show  not  only  the  relationships  among  various  variables  but  also 
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indicate  the  variability  seen  in  biological  data  of  this  kind.  These  data  are  presented 
for  use  by  those  who  wish  to  continue  development  of  mathematical  models  to  predict 
burn  injury. 


Knox,  F.  S.  HI,  Thomas  L.  Wachtel,  George  McCahan.  "Thermal  Properties  Calcu¬ 
lated  from  Measured  Water  Content  as  a  Function  of  Depth  in  Porcine  Skin" . 

In  order  to  develop  a  realistic  tissue  water  boiling  routine  for  a  mathematical 
model  of  burn  development,  it  was  necessary  to  know  the  water  content  of  skin  and 
the  thermal  properties  as  a  function  of  depth.  Split  thickness  skin  samples  were 
obtained  from  several  pigs  using  an  air  powered  dermatome.  Alternate  segments 
of  these  skin  grafts  were  processed  for  skin  water  content  determination  and  for  histo¬ 
pathologic  measurements  of  skin  thickness.  Tissue  samples  were  weighed,  dried 
and  subsequently  weighed  again  in  a  standardized  protocol  to  determine  tissue 
water  content.  In  some  instances  the  volume  of  tissue  was  also  determined  in  order 
to  allow  the  calculation  of  tissue  density.  Given  a  table  of  measured  values  of 
water  content  as  a  function  of  skin  thickness,  a  least  squares  cubic  polynomial 
was  fit  to  the  data  and  water  content  as  a  function  of  depth  was  computed  from  the 
following  formula:  w(T-d)  =  T/d  x  CW«j— +  W-p-d  where  T  is  the  total  thick¬ 
ness  of  a  slab,  W*p  is  the  fraction  of  water  computed  from  the  cubic  equation,  d  is 
the  thickness  of  the  thin  slab  at  a  depth  T-d,  and  W^-d  is  the  fraction  of  water 
above  the  thin  slab.  The  thermal  properties,  as  a  function  of  depth,  were  then 
calculated  using  the  formulations  of  Cooper  and  Trezek  (Aerospace  Med.  42(1):  24-27, 
1971.). 


Knox,  F.  S.  Ill,  Nelson  O'Young,  Daniel  D.  Reneau,  Chester  Ellis,  Thomas  L. 

Wachtel.  "BRNSIM  -  An  Analytical  Model  for  Predicting  Thermal  Cutaneous  Injury". 
BRNSIM  a  computer  model  written  in  FORTRAN  IV  run  on  a  DEC  PDP 11/40  with  80K 
of  core. 

Since,  for  all  practical  purposes,  skin  is  essentially  opaque  to  thermal  radiation 
within  a  postcrash  fire,  it  can  be  considered  to  transfer  energy  internally  by  con¬ 
duction  only.  This  energy  transfer  in  skin  can  be  described  by  the  heat  conduction 
or  Fourier  equation.  This  equation  was  implemented  in  a  12  node  model  with  appro¬ 
priate  boundary  conditions  which  were  initially:  1)  no  back  radiation  at  the  surface. 
2)  an  adiabatic  back  wall,  and  3)  no  cooling  by  blood.  The  equation  set  was  solved 
by  applying  Crank-Nicholson  method  to  the  second  order  of  partial  derivatives  and 
corresponding  explicit  methods  to  the  first  order  of  partials.  The  implicit  differ¬ 
encing  method  is  noted  for  characteristics  of  stability  and  convergence.  Correct 
increment  sizes  yield  reliable  convergence.  The  model  was  implemented  in  FORTRAN 
using  solution  techniques  of  Thomas  as  described  by  Bruce.  The  initial  model 
was  subsequently  revised  to  allow  energy  flux  across  the  surface  during  heating, 
convective  heat  loss  at  the  skin  surface  during  cooling,  heat  transfer  into  deep 
tissues  including  conduction  into  fat  and  convective  cooling  by  the  blood,  tissue 
water  boiling,  variable  thermal  properties  for  each  node  and  an  initial  temperature 
profile  from  the  surface  to  deep  sfructures.  The  model  is  run  interactively  with 
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most  variables  changeable  for  each  run.  This  model  calculates  in  depth  tissue 
temperatures  as  a  function  of  time  for  each  node.  These  temperatures  are  converted 
to  a  damage  rate  assuming  first  order  kinetics.  The  damage  rate  is  then  integrated 
to  give  an  indication  of  total  damage  at  each  node  and  the  interpolation  routine 
finds  the  depth  at  which  the  total  damage.  Cl,  equals  one.  The  parameters  of  the 
entire  model  were  adjusted  so  that  the  calculated  temperatures  agreed  with  recorded 
intraskin  temperature  profiles  and  so  that  the  calculated  damage  (threshold  depth) 
agreed  with  experimental  data  obtained  by  exposing  pigs,  as  a  human  skin  analog, 
to  simulated  postcrash  fires  of  various  intensities.  It  was  found  that  the  inclusion 
of  tissue  water  boiling,  proper  thermal  properties  for  the  skin,  an  initial  temperature 
profile  and  heat  loss  to  deep  structures  by  conduction  and  convection  were  necessary 
to  obtain  calculated  temperature  profiles  which  agree  with  the  experimentally  deter¬ 
mined  profiles.  It  was  felt  that  until  this  agreement  was  obtained,  further  adjust¬ 
ments  in  the  activation  energy  and  frequency  factors  of  the  Arrhenious  relationship 
was  an  exercise  in  futility.  The  model  prints  a  one  page  summary  of  the  calculations 
and,  at  the  user's  option,  will  produce  a  one  page  plot  of  the  calculated  temperature 
profiles  on  which  any  experimentally  determined  temperature  points  can  be  over¬ 
lay  ed.  The  model  is  designed  to  take  heat  flux  as  a  constant  value  or  as  a  table  of 
values  simulating  heat  flux  as  a  function  of  time;  thus  the  model  can  be  employed 
using  output  from  heat  flux  sensors  such  as  those  in  thermoman,  the  U.  S.  Air 
Force/ Aerotherm  instrumented  manikin. 


Knox,  F.  S.  IE,  Thomas  L.  Wachtel,  Chester  Ellis.  "A  Comparison  of  BRNSIM, 
a  Model  for  Predicting  Thermal  Cutaneous  Injury  with  Four  Earlier  Models  in  Their 
Ability  to  Precit  Exp  erim  an  tally  Determined  Burn  Injury  from  Three  Separate  Studies'''. 

BRNSIM,  a  computer  model  described  earlier,  was  run  using  activation  energies 
and  frequency  factors  as  proposed  by  Moritz  and  Henriques,  Green  and  Stoll, 

Takata,  and  Mehta  and  Wong.  Experimental  burn  data  collected  by  the  University 
of  Rochester  study  (Report  No .  394) ,  Stoll  and  USAARL  using  pigs .  humans  and 
pigs  respectively  were  used  to  test  the  various  models'  ability  to  make  accurate 
predictions  over  the  entire  range  of  burns  from  minor  epidermal  injury  to  complete 
dermal  necrosis.  Suggestions  are  made  regarding  the  use  of  such  models  and 
fabric  testing. 


Knox,  Francis  S.  in,  Chester  Ellis,  Ransom  Nock  ton .  "Thermal  Analysis  Program: 
Computer  Programs  for  Data  Base  Manipulation". 

Twelve  programs  written  primarily  in  FORTRAN  IV  are  presented  completed 
with  source  listings,  comments  and  a  brief  user's  manual.  All  of  these  programs 
have  been  run  on  a  Digital  Equipment  Corporation  PDP 11/40  with  dOK  words  of  core. 
Programs  include  the  following: 


Program  Name 


Description 


TAPE 


Reads  digitized  data  from  tape  and  constructs  named  files 
for  each  channel  of  data. 


MERGE 


SCANSEQ 


PIGBOOK 


P1GWRT 

FLUX 

TEMPFILE 

FURNACE 

PREDISC 


Reads  data  files  PIG. DAT  and  PI. DAT  on  disk  which  have 
been  previously  read  from  cards  and  merges  the  two  files 
into  the  data  base;  computes  shutter  corrections,  normalized 
burn  depths,  group  numbers,  and  converts  temperatures  to 
Centigrade . 

Scans  the  sequence  channel  for  the  tapes  read  by  program 
TAPE  and  creates  a  file  of  "events"  and  "cycles"  for  each 
tape  to  be  used  in  locating  data  for  a  given  pig. 

1)  Creates  the  pig  directory  PIG. DIR  which  locates  pig 
data  by  tape,  cycle  number,  and  burn  site. 

2)  Looks  up  data  in  the  pig  directory. 

3)  Lists  the  pig  directory  (sorted  or  unsorted)  on  the 
terminal  or  printer. 

4)  Edits  and/or  lists  the  event  tables  created  by  SCANSEQ. 

5)  Edits  and/or  lists  the  flux  directory  FLUX. DIR  created 
by  program  FLUX. 

6)  Edits  the  data  base. 

7)  Appends  the  data  to  the  data  base. 

Separates  the  digitized  data  for  a  given  channel  and  cycle 
(burn)  into  files  containing  the  data  for  a  single  event  to 
be  used  by  the  data  processing  programs. 

Computes  the  flux  seen  by  the  slug  calorimeter,  using  the 
files  created  by  PIGWRT  for  calibration  and  data  reduction. 

Computes  a  temperature  profile  from  the  copper- cons  tan  tan 
thermocouples  to  be  used  by  the  model  program. 

Computes  the  mean  furnace  wall  temperatures  seen  by  the 
Chromel-Alumel  thermocouples. 

Builds  a  file  of  values  of  variables  selected  by  input  to  be 
used  by  the  plotting  program  PLOTS  or  as  input  to  the 
discriminant  program. 
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Program  Name 

Description 

DISCRM 

Performs  a  linear  regression  analysis  on  data  from  the  data 
base  selected  by  program  PREDISC. 

PLOTS 

Plots  any  variables  in  the  PREDISC  file  vs.  any  other  on 
either  linear  or  log  scales  with  the  capability  to  parameter¬ 
ize  the  plot  by  a  third  variable  using  an  input  table  of 
values  and  symbols. 

LSTMRGD 

Builds  a  formated  file  with  labels  of  the  data  base  to  be 
printed  by  PIP  (the  system  utility  Peripheral  Interchange 
Program) . 

Knox,  F.  S.  IE,  and  Charles  Yuile.  “Pitfalls  in  the  Use  of  the  USAARL  Histopathology 
Grading  System” . 

One  hundred  biopsy  specimens  from  the  USAARL  bioassay  study  (USAARL 
Report  78-11)  were  reread  with  a  view  toward  determining  those  aspects  of  the 
grading  process  where  error  and/or  disagreement  among  the  pathologists  might 
arise.  400  biopsy  specimens  from  the  University  of  Rochester's  studies  as 
reported  in  UR  Report  338  and  553  were  prepared  and  read  to  ascertain  whether 
the  shrinkage  seen  in  more  severe  burns  in  the  USAARL  study  was  also  seen  in  the 
Univerity  of  Rochester  studies  so  that  shrinkage  correction  factors  could  be  used  to 
make  the  data  more  consistently  predictable  by  mathematical  models  such  as  USAARL 's 
BRNSIM. 


Takata,  Arthur  N.,  John  Kouss,  Thomas  Stanley,  F.  S.  Knox  III,  and  S.  C.  Knapp. 
“Thermal  Analysis  Program:  Thermal  Studies  and  Mathematical  Model  Development 
in  Support  of  USAARL  Experimental  Studies”. 

The  objective  of  this  program  is  to  develop  means  for  using  skin  simulants 
to  assess  the  thermal  protection  afforded  human  beings  in  accidental  aviation-fuel 
fires  by  various  candidate  fabrics.  In  order  to  achieve  this  objective,  USAARL 
(U.  S.  Army  Aercmedical  Research  Laboratory)  has  exposed  pigs  and  skin  simulants 
(with  and  without  fabrics)  to  aviation-fuel  fires  burning  within  a  specially  con¬ 
structed  furnace.  Results  were  then  used  by  I3TRI  (I2T  Research  Institute)  to  for¬ 
mula  te  a  criterion  of  burn  depth  and  a  code  to  translate  simulant  temperature  data 
into  burn  predictions. 

Other  supporting  endeavors  of  note  include:  a  determination  of  the  effect  of 
pig  hairs  in  shielding  the  skin  from  radiation;  a  determination  of  the  character  of 
the  beating  i.e.,  radiation  versus  convection;  a  determination  of  the  water  content 
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of  pig  tissue  with  respect  to  depth;  and  analytical  studies  indicating  the  importance 
of  thermal  properties,  flux  and  exposure  time  on  burns ,  Also  computer  codes  have 
been  prepared  to  place  all  experimental  data  on  disk  files  for  ready  computer  assess. 


Saur milch ,  Peter,  F.  S.  Knox  QI.  "User's  Manual  for  the  Burn  Prediction  Computer 
Model  Developed  by  QTRI." 

The  model  described  in  the  form  of  a  computer  program  which  takes  the  output 
of  a  physical  heat  flux  sensor  and  converts  that  measurement  of  heat  flux  to  a  pre¬ 
diction  of  skin  damage.  Included  with  this  report  are:  1)  a  description  of  model 
components,  (2)  development  of  background  mathematics,  (3)  discussion  of  the 
FORTRAN  statements  and  (4)  the  User's  Manual.  The  burn  prediction  model  is  the 
first  in  a  series  of  models  based  on  the  data  collected  using  the  USAARL  fire  simu¬ 
lation  system  and  the  porcine  bioassay  technique.  It  incorporates  tissue  water  boiling 
but  does  not  incorporate  varying  thermal  parameters,  cooling  of  the  skin  by  blood 
and  corrrection  for  thermal  shrinkage  observed  in  the  experimental  data. 


Wachtel,  Thomas  L.,  Edward  E.  Wachtel,  Francis  S.  Knox  IH,  and  Jerry  M.  Shuck. 
"The  Prevention  of  Thermal  Cutaneous  Injury  by  Hair". 

The  protective  or  detrimental  effects  of  scalp  and  facial  hair  (beards  and 
mustaches)  during  thermal  injuries  has  not  been  determined.  We  have  used  three 
separate  approaches  to  study  this  problem. 

In  a  rat  model,  half  their  flank  hair  was  alternately  clipped  or  not  clipped. 

Three  groups  of  four  animals  were  exposed  to  a  flame  heat  source  for  one,  three 
or  ten  seconds.  The  burns  were  graded  clinically  and  histologically.  Exposure 
versus  grade  curves  were  plotted.  Hair  population  and  density  calculations  were 
made.  The  protective  and  detrimental  effects  of  hair  were  determined  in  a  retro¬ 
spective  review  of  407  bum  victims  with  154  head  and  facial  bums.  A  prospective 
study  of  twenty-five  human  bums  was  undertaken  to  determine  more  accurately  the 
protective  or  detrimental  effects  of  hair  on  the  burned  skin  by  observation,  photo¬ 
graphic  record,  hair  style  and  hair  population  and  density  determinations.  The 
thermal  input  necessary  for  decomposition  and  ignition  of  hair  samples  was  determined. 
The  thermal  energy  deposited  into  the  hair  is  represented  by  the  formula  I  =  I0  exp 

x/t 

(-ox)  and  the  formula  2 •/  exp-(N*£i)»c3«  0*sin0  <10  describes  the  fraction  of  the 

incident  radiation  that  will  escape  interception  by  the  hairs.  From  the  data  and  these 
formulae  the  protective  effect  of  hair  was  determined  mathematically. 

Hair  cover  showed  significant  protection  to  the  underlying  skin  ir,  the  rat  model 
in  both  the  gross  clinical  grading  and  the  histological  sections  throughout  the  entire 
10  second  exposure  curve.  In  the  retrospective  studies  the  scalp  was  protective  in  all 
cases  except  in  the  most  severe  incincerations .  The  scalp  hair  and  facial  hair  were 
especially  protective  in  flash  exposures.  The  prospective  study  showed  beards. 
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mustaches,  sideburns  and  scalp  hair  cover  to  be  protective  of  the  underlying  skin 
in  c  .trast  to  adjacent  nude  areas  in  all  cases.  Hair  populations  and  density  deter¬ 
minations  varied  widely,  but  theoretically  and  clinically  hail  protects  the  individual 
from  burns. 


APPENDIX  C 


BSNSIM  3  -  LISTING  OF  THE  INTERACTIVE  ANALYTICAL 

MODEL  BSNSIM 


INPUTS: 


TEMPIO 

Density 

ftr 

BL 

AK 

JINC 

TEMPB 

Absorb 

Boil 

ETIME 

ITIME 

Nxtra 


Initial  Surface  Temp.  °C 
Density  of  skin  =  1.0 

Incident  Heat  Flux  either  constant  or  as  a  File  of  Fluxes 

Intemode  depth  =  22Qji 

Calculation  interval  nominally  .01  sec. 

Nodes  nominally  12 

Differences  between  TEMPIO  and  backwall  (fat/ core)  Temp.  °C 

Absorptivity  usually  .613  { 

Temperature  when  boiling  occurs  =  100.15  temperature  °C 

Exposure  Time  in  seconds 

Maximum  calculation  time  usually  80-100  sec. 

Number  of  extra  nodes  between  the  surface  and  node  1  at 
220^m 


Blood  =  Factor  to  adjust  amount  of  cooling  by  blood  usually  set  at  .001 

DEi  k  DE2  =  AE/R  from  the  Arrehemius  relationship  for  tissue  temperatures 
from  44  to  50°C  or  over  50°C  respectively 

PL]_,  PLN,  or  PLj  and  PLN2  3  log  P  =  logN  +  yloglO 

=  PL  +  PLN 

again  tor  temperatures  from  44°C-50°C  or  over  50°C 


CpCJ) 


Heat  capacity  as  a  function  of  depth 


OUTPUTS: 


IP  UTS 
TITLE 

Flux  File  I.D.,  if  used 

Flux  (I)  -  tabulated  Heat  flux  as  a  function  of  time 
DAMAGE,  W,  at  eacb  Depth  (Node) 

Maximum  Temperature 

Threshold  Depth  in  Jim 

Final  Time  -  total  calculation  time 

File  of  calculated  temperatures  for  later  plotting  by  TCPLOT 
File  summarizing  simulation 

File  of  temperature  as  printed  each  second  on  the  terminal 
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FORTRAN  IV  V01C-03D 
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MODELS, M0DEL3«M0DEL3.FTN;3 


0301 

0002 

0003 

0004 

0005 

3006 

000? 

0003 

0009 

0010 

0011 

0012 

0313 

0014 

0015 


C 

r 

C 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c. 


c. 


c 

c 


c. 

c 

c 

c 

c 

c 

c 

c 

c 

c. 


12-POINT  BURN  PREDICTION  MODEL  a**#####***#***# 

PROGRAM  PIGBRN  !  BURN  PREDICTION  MODEL  WITH  WATER  BOILING 

!  AND  USE  OF  EITHER  CONSTANT  OR  TABULATED  FLUX 
1  AND  VARIABLE  COOLING  BY  BLOOD  FROM  NODES  2  AND 
i  3  BEGINING  AT  .01  SEC  AND  LINEARLY  INCREASING 
!  TO  20  SEC  AND  THEN  REMAINGING  CONSTANT 


THIS  MODEL  WAS  DEVELOPED  UNDER  CONTRACT  FOR  THE  U.S.ARMY 
MEDICAL  RESEARCH  AND  DEVELOPMENT  COMMAND,  AND  THE  U.S.ARMY 
AEROMEDICAL  RESEARCH  LABORATORY,  FORT  RUCKER  AL.  36362, 

STANLEY  C.  KNAPP,  COL.-  MC,  COMMANDING,  BY  FRANCIS  S.  KNOX,  1 1 1,  PH.  D. 
WITH  THE  ASSISTANCE  OF  DANIEL  D.  RENEAU,  PH.D.,  NELSON  O' YOUNG, 

AND  CHET  ELLIS. M.S. 


QUESTIONS  AND  COMMENTS  SHOULD  BE  ADDRESSED  TO: 

FRANCIS  S.  KNOX,  III,  PH.D. 

ASSOC  PROF 

DEPT. ’OF  PHYSIOLOGY  AND  BIOPHYSICS 

LSU  MEDICAL  CENTER  SCHOOL  OF  MEDICINE  IN  SHREVEPORT 

BOX  33932 

SHREVEPORT,  LA.  71130 
(318)226-3134 


REAL*4  NOFIL.  ITIME 

INTEGER  BB 

LOG  I CAL* 1  RESP, YES 

’DIMENSION  T( 12) , F ( 12) , G< 12) ,H(12),Z(12),U(12),W(12) 
DIMENSION  SVC  12) ,  CP (12,  2) ,  BKC 12,  2) ,  D( 12) ,  DSCRPTC20) 
DIMENSION  ID (4) ,  FILNAMC5) ,  FLUXC 153) , SUMC 12)  -  DWC 12) 

D IMENS  ION  IFLAGC 12) ,  QC 12) ,  XTIMEC 12)  -  ZTIMEt  12) ,  JFI.AG ( 12) 
DIMENSION  WATER (12, 2), DP (12, 12) , THCONC2) , ROCONC2) , CPCONC2) 
DIMENSION  XW(3) , XTRAC8) ,  XTRALGC8) ,  XTMPC3) - XDW(8) 


EQUIVALENCE  (NOFIL, IBLNK) 

DATA  YES/'Y'/.  NOFIL/'  '  / 

DATA  MAXDIM/12/.D2/200./ 

DATA  THCON/1 .384316E-3, 4. 156843  IE-4/ 
DATA  ROC ON/6 . 1827743E-2, .33817226/ 
DATA  CPCON/.595, .405/ 


LOGICAL 

LOGICAL 

LOGICAL 

LOGICAL 

LOGICAL 

LOGICAL 

LOGICAL 


UNIT  i  INPUT  :  ' RENEAU12.DAT’ :  TABLES 

UNIT  1  SCRATCH:  ' XKRCH. TMP' :  VALUES  OF  XDW  (IF  COMPUTED) 

UNIT  2  SCRATCH:  'SCRCH.TMP';  VALUES  OF  DW  (COMPUTED) 

UNIT  3  OUTPUT  :  ' PIGBOIL.DAT' :  TEMPERATURE  PROFILES 
UNIT  4  INPUT  :  FILNAM;  NAME  OF  FLUX  FILE 

UNIT  4  OUTPUT  :  ' FIGBOIL.PLT' :  DATA  FOR  PLOTTING  TEMPERATURE 
PROFILES 

UNIT  7  OUTPUT  :  DEFAULT  'FOR007.DAT* ;  SUMMARY  PRINTOUT 
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MODEL3. MODEL3*MODEL3 . FTN  5 3 


C  ARITHMETIC  FUNCTIONS  FOR  THERMAL  CONDUCTII/ITY  (THKON)  AND 

C  DENSITY*HEAT  CAPACITY 

C 


0016 

THKON (2)  *  (THCONC 1)#Z+THC0N(2) )/<RQCON( 1)*Z+ROCON(2> ) 

0017 

r 

ROCPCZ)  -  (CPCONC 1)#Z+CPCCN(2) )/(ROCON( 1 )*Z+R0C0N(2) ) 

0018 

CALL  ASSIGN( 1, ' RENEAU12.DAT' ) 

0019 

CALL  ASSIGNC2,  'SCRCH.TMP') 

8020 

CALL  ASSIGN(3, ' PIGB0IL.DAT' ) 

c 

r 

SET  UP  PARAMETERS  FOR  THIS  RUN 

0021 

c*  *  •  « 

’  ’READ  Cl,  901) TEMP  10, DENS-  00, BL.  AK,  BOIL,  ABSORB 

0022 

901 

FORMATC7F10.5) 

8023 

REATK  1, 902) J INC, TEMPB. ITIME, ETIME,  PCWATR,  BLOOD 

0024 

902 

FORMATC 1 1 10,  7F10.5) 

0025 

READC1-901) (CP(J,2), J-l, JINC) 

0026 

READC 1, 901) CBKCJ-2), J-l, JINC) 

8027 

READ  < 1 -  90 1 ) PL2, PLN2, PL  1 , PLN 1 , DE2, DE 1 

0028 

READC 1, 901) CWATERC 1,2), 1*1,  JINC) 

0029 

CALL  CLOSE C 1 ) 

0030 

CALL  ASSIGNC 1, ' XKRCH.TMP'  ) 

0031 

FLUXCi)  -  00 

0032 

FLUXC2)  «  00 

0033 

NFLX  *  2 

0034 

FILNAMC 1)  -  NOFIL 

0035 

PPL1  -  PL1 

0036 

PPLN1  -  PLN1 

0037 

DDE1  *  DEI 

0038 

NXTRA  -  0 

0039 

NXTRA0  «  0 

0040 

700 

TYPE  910 

0041 

TYPE  911, TEMP  10, DENS, FLUXCI), BL, AK, JINC,  TEMPB,  ABSORB,  BOIL 

0042 

TYPE  104, PL2, PLN2, PL  1 , PLN 1 , DE2, DE 1 , ETIME, ITIME, NXTRA, BLOOD 

0043 

IF  C NXTRA)  TYPE  7003, CXTRAC I) , 1-1, NXTRA) 

0045 

7003 

FORMATC '  0EXTRA  NODES:  ',8F6.D 

0046 

4 

TYPE  906 

0047 

906 

FORMATC IX, 'CONTINUE  ?  Y/N',T15,$) 

0048 

ACCEPT  220, RESP 

0049 

IFCRESP.NE.YES)  GO  TO  290 

0051 

5 

CONTINUE 

0052 

TYPE  699 

8053 

693 

FORMATC'  1  FUNCTION  *"S'/'81  -  REENTER  INITIAL  YALUES'  / 

1, '  2  -  CHANGE  SELECTED  INITIAL  VALUES'/'  3  -  NO  CHANGES' / 

2, '  4  -  EXIT' //' SENTER  FUNCTION  NO.:') 

0354 

ACCEPT  801-ANSR 

0055 

IANSR-ANSR 

0056 

TYPE  710  , ANSR 

3057 

710 

FORMATC ' 0ANSR*'  ,  F 10 . 0) 

vO 

BB  -  0 

8059 

IF ( IANSR . EQ.2) BB  *  1 

0061 

GO  TO  (20,80,230.  lill) IANSR 

0062 

20 

TYPE  800, TEMP  10 

0063 

800 

FORMATC//"  PLEASE  PROVIDE  THE  REQUESTED  I'ALUES:'// 

1T5-  '  TEMP  10'  ,  G16.3.  T15,  $) 

0064 

ACCEPT  801, TEMP  10 

0065 

801 

FORMATCF10.0) 

0066 

GO  TO  (200) BB 

0067 

21 

TYPE  802- DENS 

86 
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0068 

802 

F0RMATC/T5, 'DENS' .  G16.8,  T15,  $) 

0069 

ACCEPT  801/ DENS 

0070 

. 

GO  TO  (200) BB 

0071 

22 

TYPE  1 

0072 

1 

FORMAT (✓' fFLUX  FILE  NAME  (BLANK  IF  CONSTANT  FLUX):  ') 

0073 

ACCEPT  2/FILNAM 

0074 

2 

FORMAT (5A4) 

0075 

IF  (FILNAMC 1) -EQ.NCFIL)  GO  TO  799 

0077 

CALL  ASSIGN(4/ FILNAM) 

0078 

READ (4) ID/TDELT/NFLX 

*079 

TYPE  3/  ID/TDELT/NFLX 

n 

UUuu 

READ (4) (FLUXCI),  I-l.NFLX) 

0081 

CALL  CLOSE  (4) 

0032 

GO  TO  (200) BB 

0083 

GO  TO  23 

0084 

799 

TYPE  788/FLUXd) 

0085 

788 

FORMATS  ^CONSTANT  G-l/ALUE  -  d  616.8. ‘NOH.  IT  SHOULD  BE  Q 

0086 

ACCEPT  80 1 . FLUX( 1 ) 

0087 

FLUX(2)«FLUX(1) 

0088 

NFLX  -  2 

0089 

DO  777  1-1.4 

0090 

777 

ID  ( I )'  -  IBLNK 

8091 

GO  TO  (200) BB 

0092 

23 

TYPE  804. BL 

0093 

804 

F0RMAT(/T5, ' BL' . G16.8, T15. S) 

0094 

ACCEPT  80 1 . BL 

0095 

GO  TO  (200) BB 

0096 

24 

TYPE  805, AK 

0097 

805 

FORMAT (/T5, 1 AK' , G 16.8- TI5,  *) 

0098 

ACCEPT  801. AK 

0099 

GO  TO  (200) BB 

0100 

25 

TYPE  806- J INC 

0101 

806 

FORMAT(/T5, ' JINC'  . I10.T15.S) 

0102 

ACCEPT  820. J INC 

0103 

IF  (JINC.GT.MAXDIM)  JINC-MAXDIM 

0105 

GO  TO  (200) BB 

0186 

820 

FORMAT( I 10) 

0107 

26 

TYPE  807/TEMPB 

8108 

807 

FORMAT  (/T5, 'TEMPB' - F10.0,  T15,  $) 

0109 

ACCEPT  801 -TEMPB 

0110 

GO  TO  (200) BB 

0111 

27 

TYPE  808.ETIME 

0112 

808 

FORMAT! /T5, ' ETIME' . G 16.8,  T15.  $) 

0113 

ACCEPT  801 .ETIME 

0114 

GO  TO  (200) BB 

0115 

28 

TYPE  809, PL  1 

0116 

809 

FORMAT! /T5, '  PL l '  ,  G16. 8.  T15  -  $) 

0117 

ACCEPT  801.PL1 

0118 

PPL1-PL1 

0119 

GO  TO  (200) BB 

0120 

29 

TYPE8 10- PLN 1 

0121 

810 

FORMATUT5,  'PLN1'  -  G16.8,  T15,  $) 

0122 

ACCEPT  801, PLN 1 

0123 

PPLN1-PLN1 

0124 

GO  TO  (200) BB 

0125 

38 

TYPE  811, PL2 

0126 

811 

FORMAT ( /T5 - '  PL2 ' ,  G 1 6 . 8 ,  T 1 5 ,  $) 

0127 

ACCEPT  801.PL2 

) 
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0128  GO  TO  (200) BB 

0129  31  TYPE  812.PLN2 

0130  812  FORMATC/T5. 'PLN2'  ,G16.8.T15.$) 

0131  ACCEPT  801.PLH2 

0132  GO  TO  (208) BB 

0133  32  TYPE  8 13. DEI 

0134  813  F0RMAT(/T5» '  DEI'  ,G16.8.T15,$) 

0135  ACCEPT  801, DEI 

0136  DDE  1 -DEI 

0137  GO  TO  (200) BB 

0138  33  TYPE  814, DE2 

0139  814  F0RMATC/T5, ' DE2'  ,G16.8.T15.S) 

0140  ACCEPT  801, DE2 

0141  GO  TO  (200) BB 

0142  34  TYPE  815. ITIME 

0143  815  FORMAT (/T5, ' ITIME'.G16.8,T15,$) 

0144  ACCEPT  801, ITIME 

0145  GO  TO  (200) BB 

0146  35  TYPE  816, ABSORB 

0147  816  FORMAT(/T5, 'ABSORB', G16. 8,  T15.S) 

0148  ACCEPT  801, ABSORB 

0149  GO  TO  (200) BB 

0150  36  TYPE  817. BOIL 

0151  817  FORMAT(/T5, 'BOIL' . G16.8,  T15.  $) 

0152  ACCEPT  801, BOIL 

0153  GO  TO  (200)  BB 

0154  7100  TYPE  7005. NXTRA 
0155  NXTRA0  ■  NXTRA 

0156  7005  F0RMAT(/T5. 'NO.  OF  EXTRA  NODES  ',I4,2X,$) 

0157  ACCEPT  7007, NXTRA 

0158  7007  FORMAT( 14) 

0159  IF  (NXTRA. EQ.0)  GO  TO  7230 

0161  IF  (NXTRA. GT. 8)  NXTRA -8 

0163  IF  (NXTRA0)  TYPE  7009, (XTRA( I) ,  1-1, NXTRA0) 

0165  7009  FORMATS  ©CURRENT:  '.8F6.D 

0166  TYPE  7011 

0167  7011  FORMAT ( '  0ENTER  NEW  (VALUES  SEPARATED  BY  COMMAS,  OR  <CR>  IF  PROGRAM' 

1/'  IS  TO  CALCULATE  WLUES:') 

0168  ACCEPT  70 13. XTRA 

0169  7013  FORMAT(8F8. 1) 

0170  IF  (XTRA(D)  GO  TO  7210 

C 

C  NUMERATOR  IN  NEXT  STATEMENT  IS  SPECIFIC  FOR  N-POINT  MODEL 
C 

0172  DXTRA  -  D2/(NXTRA+1) 

0173  DO  7200  I- 1, NXTRA 

0174  7200  XTRA ( I )  •  DXTRA* I 

0175  7210  DO  7220  I «1, NXTRA 

.,76  7220  XTRALG(I)  -  ALOG(XTRA( I) ) 

0177  7230  CONTINUE 

0178  NXTRA0  -  NXTRA 

0179  GO  TO  (200) BB 

0180  37  TYPE  8 18. BLOOD 

0181  818  FORMATS' $  BLOOD  ■  '.F10.4.3X) 

0182  ACCEPT  801, BLOOD 

0183  GO  TO  200 

0184  80  TYPE  600 

0135  600  FORMAT(/T5, 'PICK  A  NUMBER',  /VT10, ' 1 -TEMP  10', T30. ' 8-ETIME' . // 
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1T10, * 2-DENS' - T30, ' 9-PL1' . /T10, ' 3-Q1' , T30, * 10-PLN1',// 

1T10. ' 4-BL'  -  T30, ' 1 1 -PL2' -  /TT0,  ' 5-AK'  -  T30. '  12-PLN2* . // 

1T10, ' 6-J INC' , T30, ' 13-DE1'  - /T10. ' 7-TEMPB' , T30.  ' 14-DE2' ,  ✓✓ 
1T10, '  15-  ITIME'  , T30, '  16-ABSORBTIWITY'  / 

1T10. ' 17-BOIL' .  T30. *  18-EXTRA  NODES'  /T10. ' 19-BLOOD' ,  T30/V) 

0186 

BB-1 

0187 

ACCEPT  601, INUM 

0188 

601 

FORMAT(I2) 

0189 

r 

GO  TO  (20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35.36 
1,7100,37) INUM 

c 

ANOTHER  WALUE? 

0190 

Q 

200 

TYPE  210 

0191 

210 

FORMAT (/'SANOTHER  CHANGE  ?  TrTE  Y/N' ) 

0192 

ACCEPT  220, RESP 

0193 

220 

FORMAT(Al) 

0194 

IF  (RESP.  EG).  YES)  GO  TO  80 

0196 

GO  TO  700 

0197 

230 

REWIND  1 

0198 

TYPE  907 

0199 

907 

FORMAT( ' 0ENTER  MODEL  NAME  OR  DESCRIPTION  (80  CHARACTERS)') 

0200 

ACCEPT  908, DSCRPT 

0201 

908 

FORMAT(20A4) 

0202 

WRITE (3.  909) DSCRPT 

0203 

WRITE (7.  909) DSCRPT 

0204 

909 

FORMATS  MODEL  NAME  OR  DESCRIPTION:  ' , 20A4) 

C 

C  ASSIGN  LOGICAL  UNIT  4  FOR  PLOT  OUTPUT 


0205 

C 

CALL  ASSIGN(4.  '  PIGBOIL.PLT' ) 

0206 

WRITE (4) DSCRPT 

0207 

AJ-JINC 

0208 

Q1  -  FLUX(l) 

0209 

r 

H 1 «BL/(AJ- 1.0) 

L 

c 

INITIALIZE  DEPTH  NODES  D(J) 

0210 

D(l)  -  -16. 

0211 

DO  106  I -2, JINC 

0212 

D( I)  -  H1*(I-1)*1.E4 

0213 

106 

DC  I )  »  ALOG(D<I)) 

3214 

DTJ  -  TEMPB/(JINC-1) 

0215 

D06J-1.  JINC 

0216 

WATER (J. 1)  -  WATER(J, 2) 

0217 

CPCJ. 1)  -  CP( J, 2) 

0218 

BK(J.  1)  «  BK( J,  2) 

0219 

XTIME(J)  -0. 

0220 

ZTIMECJ)  -0. 

0221 

IFLAG(J)  -0. 

0222 

JFLAG(J)  -0. 

0223 

6 

T(J)  -  DTJ*(J-1)+TEMPI0 

0224 

WRITE (4)0. , (T( J) . J-l, J INC) 

0225 

WRITE(3. 910) 

0226 

WRITEC7,  910) 

0227 

910 

FORMAT ( 1H0, 55X. 'SKIN  DIFFUSION  DATA' A 55X. ' INPUT  PARAMETER 

0223 

WRITE (3,  91 1 ) TEMP  10, DENS, Ql, BL. AK, J INC.  TEMPB- ABSORB.  BOIL 

0229 

WR I TE ( 7 .  9 IDTEMPI0. DENS.  Ql.BL.AK,  JINC. TEMPB.  ABSORB.  BOIL 

89 
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0230 


0231 

0232 

0234 

0235 

0237 


0238 

0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0258 


0259 

0261 

0262 

0263 

0264 

0266 

0267 


0269 

0271 

0272 

0273 


911  FORMAT ( IX. 52X.  8H  TEMPI0-. G14.6/54X. 5HDENS-. G14.6/S4X,  3HQ1 «. G14. 6/ 
154X. 3HBL-,  G 14. 6/54X.  3HAK-, G 14. 6/54X.  5HJ INC-. 16/S4X. 6HTEMPB-.  F8 . 4/ 
154X.  'ABSORB-'  .G14.S/54X.  'BOIL-' .G14.6/) 

WRITE(3. 104) PL2. PLN2. PLl. PLN1,  DE2. DEI. ETIME.  ITIME. NXTRA.  BLOOD 
IF  (NXTRA)  WRITEC3.7003) (XTRA(I).  I-l. NXTRA) 

WR ITE (7.  104) PL2.  PLN2. PL  1 .  PLN l , DE2. DE 1 . ETIME. ITIME. NXTRA. BLOOD 
IF  (NXTRA)  MR ITE (?. 7003) (XTRA( I) , I-l. NXTRA) 

104  FORMAT! IX.  53X. '  PL2-'  . F10.2/S4X. ' PLN2»' » F10.2/54X. ' PL1-' . F10.2/ 
154X.  'PLN1-' .F10.2/54X. 'DE2-'  -F10.2/54X, ' DE1-' . F10.2/ 

154X. 'ETIME-', F 10. 2/S4X.' ITIME-' , F10.2/S4X, *  NXTRA-' , I4/S4X 
1,  'BLOOD-'  .F 10. 4) 

WRITE(3, 3) ID. TDELT, NFLX,  (I.FLUX(I),  I-l.NFLX) 

WRITE(7.  3) ID, TDELT,  NFLX, ( I,  FLUX( I) ,  I-l.NFLX) 

3  FORMAT (' 0FLUX  FILE  I.D.:  ' , 4A2, F7 .2. 14://'  FLUX(I)-' 

l/('  ',  10(15. F8. 3))) 

JJJJ-0 

F(l) — BK(2. 1)/(2.0*H1*HI)-BK( 1, 1)/(2.0*H1*H1) 

G( 1)  -  (BK(1. l)+BK(2, 1))/(2.0*H1*H1)+DENS*CP(1. 1)/AK 
H( 1) -0.0 
ITTR  -  0 
IFLX  -  1 

EITIM1  -  ITIME+1 . 

IF  (FILNAM(l) .EQ.NOFIL)  TDELT-AK 
FFDG  -  TDELT/AK 
KFDG  -  FFDG+.0001 
TMPMAX  -  0. 

QCONST  -  ABSORB*60.892 
BLUD  -  0. 

M  -  -1 
TIME  -  0. 

ITTLG  -  0 
GO  TO  12 
C 

C  STATEMENT  11  AUTOMATICALLY  CHOOSES  PROPER  INTERIM  IN  FLUX  TABLE 
C  FOR  COMPUTATION  OF  QT  AND  Ql  FOR  EITHER  CONSTANT  OR  TRIABLE 
C  FLUX. 

C 

C  KFDG  (-FFDG)  •  1  FOR  CONSTANT  FLUX 

C  -  C RATIO)  OF  TABULAR  TIME  STEP  TO  MODEL  TIME  STEP  FOR 

C  TABULATED  FLUX 

C  C...]  MEANS  INTEGER  WLUE  OF _ 

C 

11  IF  (MOD( ITTR. KFDG) .EQ.0. AND. IFLX. LT. NFLX)  IFLX-IFLX+t 
ITTR  -  ITTR+1 

P  -  ( ITTR-KFDG#( IFLX-2) ) /FFDG 

QT  -  ( 1 . -P) *FLUX( IFLX-1 ) +P*FLUX( IFLX) 

Ql  •  QT*QCONST 
JJJJ  •  JJJJ+1 
T!ME-JJJJ*AK 

IF  (TIME. GE.. 01. AND. TIME. LE. 20.)  BLUD- (TIME-,0 t)/19.99*BLOOD 
C 

C  WLUE  23.9  IS  QUESTIONABLE  SHOULD  BE  STUDIED  (RESOLVED) 

C 

IF(TIME.GE. ETIME)  Ql— 5.E-4»'(7(  D-23.9) 

Z(D— F(1)*T(2)-((BK(1.  1)+BK(2,  15 )/(2.0*Hl*Hl)-(DENS*CP( 1.  1))/AK) 
l*T(l)+Ql 
N-J INC- l 
DO10J-2.N 
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0274 

0275 

0276 

0277 

0278 

0280 

0281 

0282 

0283 

0284 

0286 

0287 

0288 

0289 

0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 

0298 

0300 

0301 

0303 

0305 

0306 

0308 

0309 

0310 

0311 

0312 

0313 

0314 

0315 

0317 

0318 

0319 

0323 

0321 

0322 

0323 


0324 

0326 


0328 

0329 

0330 

0331 


F(J)—  BKCJ+l,  1)/(2.0*H1*H1) 

GCJ)“(BK(J,  1)+BK(J+1.  1))/'C2.0*H1*H1)+DENS*CPCJ,  1)/AK 
H(J)—  BKCJ.  1)/<2.0*H1#H1) 

ZCJ) “-F(Jj*TCJ+l)-( (BKCJ,  D+BKCJ+l,  1))/(2.*H1*H1)-DENS*CP(J,  1) 

i/ak)*t(J)-h(J)*tcj-i) 

IF  (J  GT.3)  GO  TO  10 

ZCJ)  -  ZCJ)-t.67S#Hl/'BK(J,  1)*BLUD*(TCJ)-TEMPI0+TEMPB) 

10  CONTINUE 

FCJINC) "0.0 

GCJ INC) “(BKCJ INC.  D+BKCJINC-1,  1)  )/'(2.0*Hi*Hl)+DENS#CP(JINC.  1)/AK 
H(J INC) “-(BKCJ INC, 1)+BK(JINC-1. 1) )/(2.0*Hl*Ht) 

DT“T(JINC-1)-(TEMP I0+TEMPB) 

Z( JINC) «(H( J INC) +(DENS#CP ( JINC, 1)/AK) )*TCJINC) 

1-H ( J INC) *T( J INC- 1) -BKCJ INC. 1)*DT/H1**2 
WC 1) “G( 1) 

U(1)-2C1)/W(1) 

DO40J-2,  JINC 
JMl-J-l 

SWCJMl)  “FCJMD/'WCJMD 
W(J)-GCJ)-H(J)*S^(JM1) 

40  UCJ)-(2(J)-H(JXI(JM1))/W(J) 

TCJINC) “UCJINC) 

KK-JINC-1 

DO50J-1.KK 

KMJ-JINC-J 

IF  ( IFLAGCKMJ) .EQ. 1)  GO  TO  203 
TCKMJ)  “U  (KMJ) -SWCKMJ) *T(KMJ-H ) 

IF  (JFLAGCKMJ) .EQ. 1)  GO  TO  50 
IF  CTCKMJ) .LT.BOIL)  GO  TO  50 
TCKMJ)  -  BOIL 
IF  (KMJ.EQ. 1)  GO  TO  109 
Q(KMJ)  “  BKCKMJ,  1 ) *(TCKMJ) -TCKMJ+1 ) )/Hl 
GO  TO  205 
109  QCKMJ)  -  QT 

205  XTIMECKMJ)  -  539. *H1 /QCKMJ )*WATERCKMJ. 1) 

ZTIMECKMJ)  -  XTIMECKMJ) +TIME 
IFLAGCKMJ)  -  1 
GO  TO  50 

203  IF  (TIME. LT. ZTIMECKMJ) )  GO  TO  50 
WATER (KMJ. 1)  ■  PCWATR 
CPCKMJ,  1)  =  ROCPC  WATER  (KMJ.  D) 

BKCKMJ,  1)  •  THKONC  WATER  (KMJ,  D) 

IFLAGCKMJ)  -  0 
XTIMECKMJ)  -  0. 

JFLAGCKMJ)  -  1 
50  CONTINUE 
C 

C  INTERPOLATE  EXTRA  TEMPERATURES  BETWEEN  SURFACE  AND  END  NODE  IF  CALLED  FOR 
C 

IF  (NXTRA.EQ.0)  GO  TO  7360 
IF  (TC2) .NE.TC l) )  GO  TO  7320 
C 

C  CONSTANT  TEMPERATURE 
C 

DO  7310  I ■ 1 , NXTRA 
7310  XTMPCI)  -  TC2) 

GO  TO  7360 

7320  IF  (TC2).NE.TC3))GO  TO  7340 
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C 

C  LINEAR  INTERPOLATION 
C 


0333 

DO  7330  I ■ l , NXTRA 

0334 

P  -  XTRACn/DE  IDCl)  •  0. 

0335 

7330 

XTMPCI)  -  ( 1 .-P)*T< 1)+P*TC2) 

0336 

r 

GO  TO  7360 

C  3-POINT  LAGRANGE  INTERPOLATION  FOR  EQUALLV  SPACED  ABSCISSAE 

0337 

7340 

DO  7350  I -1, NXTRA 

0333 

P  •  CXTRA(I)-D2)/D2  !D(1)  «  0.  (SURFACE) 

0333 

7350 

XTMPCI)  •  ,5*P*(P-1 . )*T( l)+(  1  .**P**2)*T(2)  +.5*P#(P+1 . )*T(3) 

0340 

7360 

CONTINUE 

0341 

IF  (ABSCETIME-TIME) .GT..5*AK)  GO  TO  48 

0343 

DO  44  I-l.JINC 

0344 

IF  (IFLAG(I) .£0.3)  GO  TO  44 

0346 

WATER C I. 1)  -  (ZTIME(I)-TIME)/XTIME(I)*(WATER(I. 1)~PCWATR)+PCWATR 

0347 

CPCI.  1)  -  ROCPCWATERC I.  D) 

0346 

BK(  I.  1)  -  THKONC WATER ( I.  D) 

0349 

44 

CONTINUE 

0350 

DO  45  I-l.JINC 

0351 

XTIMECI)  •  0. 

0352 

IFLAG(I)  -  0 

0353 

45 

JFLAG(I)  -I 

0354 

48 

CONTINUE 

0355 

IF  (TCI) .GT.TMPMAX)  TMPMAX-TC1) 

0357 

ITFLG  -  -1  ! ITFLG  SET  TO  0  IF  ANV  TEMPERATURE  .GE.  44  DEGREES 

0358 

DO  13  J-l.JINC 

0359 

IF  (T(J) .GE.44. )  GO  TO  14 

0361 

DWCJ)  -  0. 

0362 

GO  TO  13 

0363 

14 

CONTINUE 

0364 

ITFLG  -  0 

0365 

IF ( T( J ) . GE.50. )  GO  TO  51 

0367 

PH  •  PPL  1 

0368 

PLN 1  -  PPLNl 

0369 

DEI  -  DDE1 

0370 

GO  TO  15 

0371 

51 

PL1-PL2 

0372 

PLN1-PLN2 

0373 

DE1-DE2 

0374 

15 

DWLH-PLl+?LNl-DEl/(TCJ)+273.) 

0375 

DW( J)  -  EXPCDWLN) 

0376 

13 

CONTINUE 

0377 

WRITE(2)DW 

0378 

IF  (NXTRA. EQ.0)  GO  TO  7400 

0730 

DO  7395  J-l. NXTRA 

J381 

IF  (XTMP(J) .GE.44.)  GO  TO  7370 

0333 

XDW(J)  •  0. 

0384 

GO  TO  7395 

0385 

7370 

IF  (XTMP(J) .GE.50.)  GO  TO  7380 

.0387 

PL1  «  PPL  1 

0388 

PIN!  -  PPLNl 

0389 

DEI  •  DDE  1 

0390 

GO  TO  7390 

0391 

7380 

PL!  •  P12 

0392 

PLN  t  •  PLN2 
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0393 

0394 

0395 

0396 

0397 

0398 

0399 

0480 

0402 

0404 

1,-406 

0407 

0408 

0409 

0410 

0411 

0412 

0413 

0415 

0416 

0417 

0418 

0419 

0420 

0421 

0423 

0424 

0425 

0426 

0427 

0428 

0429 

0430 

0431 

0432 

0433 

0434 

0435 

0436 

0437 

0439 

0440 

0441 

0442 

8443 

6444 

0445 

0446 

8447 

0440 

0449 

0450 

0451 


0452 

0453 


DEI  •  DE2 

7390  DWLN  •  PLl+PLNl-DEl/(XTMP(J)+273.) 

XDW(J)  ■  EXP (DWLN) 

7395  CONTINUE 

WRITEmXDU 
741 0  CONTINUE 

EMTIME  -  AINTC 1000. *(TIME+. 00301) )/100. 

IF  (TIME. IT. 10. . AND. AMOD( EMTIME.  1.) .EQ.0. .OR. TIME. GE. 10. .AN£. 

1  AMOD (EMTIME, 10.) .EQ.0.)  WRITE(4)TIME. T 
IF  (ITFLG. AND. TIME. GE.ETIME)  GO  TO  12 
IF(JJJJ.EQ.M*100.OR.JJJJ.EQ.l)  GO  TO  12 
GO  TO  11 

12  WRITE (3. 840) TIME 

TYPE  840. TIME.  (T( I) . XTIMEC I) .  I-l . J INC) 

840  FORMATS  0'  ,  45X.  5HTIMEa,  F18.6: .  T4.  *T-  * .  6X. '  XTIME-'  /V  '.2612.4)) 
WRITE (3. 842) CXTIME(I). I-l.JINC) 

842  F0RMAT(2X. '  XTIME"'  . F10.5) 

WRITF.(3. 9t4)T(l).  CP(  1.  1).BK(1.  1) 

IF  (NXTRA.EQ.0)  GO  TO  7428 
DO  7410  J-l.NXTRA 
7410  WRITE (3. 914)XTMP(J) 

7420  WRITEC3.914) (T(J).CP(J. 1).BK(J. 1). J-2. JINC) 

914  F0RMAT(2X. ' T-* . G16.5: >2X. ' CP*' . G16.5. 2X. ' BK-' , G16.5) 

M-M+l 
EM  ■  M 

IFCTIME.GE. ITIME, OR. ITFLG. AND. TIME. GE.ETIME)  GOTO  100 
GO  TO  11 
100  CONTINUE 
REWIND  2 
READ(2)DW 
DO  17  I  «*  1 ,  J  INC 
17  SUM(I)  ■  .5*DW(I) 

DO  16  J-2.JJJJ-1 
READ(2)D1-! 

DO  16  I-l.JINC 
16  SUM( I ) »SUM( I ) ♦DW( l ) 

READ(2)DW 
DO  88  I-l.JINC 

88  W(I)  -  (SUM( I)+.5*DW( I) )*AK 
REWIND  2 

IF  (NXTRA.EQ.0)  GO. TO  7460 
REWIND  1 
READ ( 1 ) XDW 
DO  7430  J-l.NXTRA 
7430  SUM(J)  -  .5*DW(J) 

DO  7440  I-2.JJJJ-1 
READ ( 1) XDW 
DO  7440  J-l.NXTRA 
7440  SUM(J)  -  SUM( J) +XDW( J) 

READ(1)XDW 
DO  7450  J-l.NXTRA 

7450  XW(J)  -  (SUM(J)*.5*XDW(J)1#AK 
REWIND  1 
7460  CONTINUE 
C 

C  SELECT  U(J)  AND  D(J)  NEAR  U(J)  -1 
NN  •  3 

DO  110  J-l.JINC 
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0454  JLT1  -  J 

0435  IFCW(J) .GT.l.)GO  TO  110 

0457  IF(W(J).EQ. l.)GO  TO  111 

0459  IF  CJ.EQ.l)  JLT1-2 

0461  I?  CJ.EQ.JINC)  JLT1-JINC-1 

0463  GO  TO  112 

0464  110  CONTINUE 

0465  JLT1  *  JINC  !  EXTRAPOLATE 

0466  NN  •  2 

0467  112  WRITEC3.  120HWXK) .  K-JLT1-1.  JLT1+1) 

0468  WRITEC?,  120) CWCK) ,  K-JLT1-1.  JLT1+1) 

0469  120  FORMATC/C1X.  *W«' .E20.5)) 

0470  WRITEC3. 121) CDCK) .K-JLT1-1. JLT1+1) 

0471  WRITEC7.  121) CDCK).K*JLTI-1.  JLTl+i) 

0472  121  FORI’flTC/'C  IX. '  D*'  ,  E20.5) ) 

0473  IF  (NXTRA0.EQ.0.OR.JLTI .GT.2)  GO  TO  1005 

0475  TYPE  70 IS 

0476  7013  FORMATC'  0W*1  LIES  ABOWE  NODE  2  .  INTERCOLLATING  l/ALUES  OF  D  AND  W  / 

1.'  COMPUTED  FROM  INTERPOLATED  WLUES  OF  D  AND  TEMPERATURE'/') 

0477  WRITEC3.7015) 

0478  WRITEC?. 7015)  ' 

0479  WRITEd)D(l) 

0480  WRITEC2)W<1) 

0481  DO  7470  J-l.NXTRA 

0482  WRITEC DXTRALGCJ) 

0483  7470  WRIT£C2)XWCJ) 

0484  DO  7480  J-2.JINC 

0485  WRITEC 1) D(J) 

0486  7480  WRITEC2)W(J) 

0487  REWIND  1 

0488  REWIND  2 . 

0489  DO  7490  J-l.JINC 

0490  READ ( 1 ) D ( J ) 

0491  7496  READ(2)W(J) 

0492  REWIND  1 

0493  REWIND  2 

0494  NXT3A0  •  0 

0495  GO  TO  7460 

0496  1005  CONTINUE 

0497  NXTRA0  •  NXTRA 

0498  IF  CWCJLTH-l)  .EQ.0.  .AND.NN.EQ.3)-  NN-2 

0580  CALL  DEPTHCDC JLTl- l).W(JLTl-l). NN. TD. IERR) 

0531  IF  (IERR.EQ.0)  GO  TO  102 

0503  TYPE  1002 

0504  WRITEC3.  1002) 

0505  WRITEC7. 1002) 

0506  1002  FORMAT!* OERROR  IN  SUBROUTINE  ’DEPTH*.  EXITING.'/*) 

*507  GO  TO  290 

0509  102  CONTINUE 

0509  IF  CNN.EQ.2.AND.JLT1.EQ. JINC)  GO  TO  1000 

0511  GO  TO  1020 

0512  til  TD*EXPcDcJ>) 

0513  999  FORMAT! A  IX.  ' THRESHOLD  DEPTH  *  *.620.4) 

0514  GO  TO  1020 

0515  1000  WRITEC3.  100DMAXDIM 

0516  WRITEC?.  IQODMAXDIM 

351?  100 1  FORMATC/*lX, ‘THE  MODEL  BLEW  UP:  DAMAGE  >  i  AT  NODE  * . I2/> 

0518  1020  W8ITEC3.  lOl)CU(I). 1*1. JINC) 
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0519 

0520 

101 

0521 

2010 

0522 

0523 

0524 

0525 

0526 

0527 

0528 

1021 

0529 

103 

0531 

0532 

C 

0533 

C 

0534 

290 

0535 

300 

0536 

0537 

0538 

98 

0539 

0540 

0542 

mi 

0543 

0544 

0545 

0546 

054S 

1112 

0549 

0550 

99 

0551 

0552 

tO50 

0553 

OS  4 

1051 

0555 

0556 

WRITEC7, 2010)  CWCD.EXPCDCD).  1-1,  JINC) 

FORMATC/UX.  *W-*  ,E20.5>) 

FORMATC/C IX, ' W  «' . E20.S.  5X. ' AT  DEPTH  (IN  MICRONS)-' ,G20. 6) ) 

.  WKITEC3, 103)TMPMAX 
WRITE (7,  103) TMPMAX 
WRITE<3.999)TD 
WR ITE(7 , 999) TD 
WRITEC3.  1021) TIME 
WRITE  (7,  102DTIME 
FORMATCA  IX, 'FINAL  TIME  -  *,F7.2) 

FORMAT C'0MAX I MUM  TEMPERATURE  •  '.F8.3) 

TYPE  999, TD 
TYPE  103. TMPMAX 
TYPE  1021.  TIME 

CALL  CLOSE  (4) 

TYPE300 

FORMATC/' SDO  YOU  WANT  ANOTHER  RUN  7  TYPE  Y/N' ) 

WRITEC3. 98) 

WR1TEC7,  98) 

FORMAT  ('  1'  ) 

ACCEPT  220. RESP 
IF  ( RESP. EQ. YES)  GO  TO  700 
CONTINUE 
IOER  •  0 

CALL  DELEFT ( 1,  IOER) 

CALL  DELEETC2. IOER) 

IF  (IOER)  TYPE  1112 

FORMATC ' 0ERROR  IN  DELETING  FILECS)  "SCRCH.TMP"  OR' 

I.'  "XKRCH.TMP"  '  ) 

TYPE  99 

FORMATC' 1' .T20. 'PRINTED  OUTPUT  IN  FILE  "  PIGBOIL.DAT' "/'T29. 

1' USE  *PIP*  TO  PRINT'  ) 

WRITEC5. 1050) 

FORMATC ' 0' .T20.' TEMPERATURES  FOR  PLOT  IN  "  PIGBOIL. PUT  "  /T29. 
I'  USE  TCPLOT  AND  THEN  C7.  U 3SOSM  TO  PLOT') 

TYPE  1051 

FORMATC ' 3' . T20.  ‘  SUmpp  •  PRINTOUT  IN  FILE  FORQOF.DAT”  / 
r  ' , T27. ' USE  PIP  TO  PRINT'/) 

CALL  EXIT 
END 
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FORTRAN 

IV 

STORAGE  MAP 

NAME 

OFFSET 

ATTRIBUTES 

T 

000006 

REAL*4 

ARRAV  (125 

F 

000066 

REALM 

ARRAY  (12) 

G 

000146 

REALM 

ARRAY  (12) 

H 

000226 

REALM 

ARRAY  (12) 

Z 

000306 

REALM 

ARRAY  (12) 

U 

000366- 

REALM 

ARRAY  (12) 

U 

000446 

REALM 

ARRAY  (12) 

SV 

000526 

REALM 

ARRAY  (12) 

CP 

000606 

REALM 

ARRAY  (12,2) 

VECTORED 

BK 

000746 

REALM 

ARRAY  (12,2) 

VECTORED 

D 

001106 

REALM 

ARRAY  (12) 

DSCRPT 

001166 

REALM 

ARRAY  (20) 

ID 

001306 

INTEGERS 

ARRAY  (4) 

FILNAM 

001316 

REALM 

ARRAY  (5) 

FLUX 

001342 

REALM 

ARRAY  (150) 

SUM 

002472 

REALM 

ARRAY  (12) 

DU 

002552 

REALM 

ARRAY  (12) 

I  FLAG 

002632 

INTEGERS 

ARRAY  (12) 

Q 

002662 

REALM 

ARRAY  (12) 

XTIME 

002742 

REALM 

ARRAY  (12) 

ZTIME 

003022 

REALM 

ARRAY  (12) 

JFLAG 

003102 

INTEGERS 

ARRAY  (12) 

WATER 

003132 

REALM 

ARRAY  (12,2) 

VECTORED 

DP 

003272 

REALM 

ARRAV  (12- 12) 

VECTORED 

THCON 

004372 

REALM 

ARRAY  (2) 

ROCON 

604402 

REALM 

ARRAY  (2) 

CPCON 

004412 

REALM 

ARRAY  (2) 

XW 

004422 

REALM 

ARRAY  (8) 

XTRA 

004462 

REALM 

ARRAY  (8) 

XTRALG 

004522 

REALM 

ARRAY  (8) 

XTMP 

004562 

REALM 

ARRAY  (8) 

XDW 

004622 

REALM 

ARRAY  (8) 

• 

NOFIL 

004662 

REALM 

VARIABLE 

ITIME 

011534 

REALM 

TRIABLE 

BB 

011540 

INTEGERS 

TRIABLE 

RE5P 

011542 

LOG ICAL*1 

VARIABLE 

VES 

004666 

logical*i 

TRIABLE 

IBLNK 

004662 

INTEGERS 

VARIABLE 

MAXDIM 

004670 

INTEGERS 

VARIABLE 

D2 

004672 

REALM 

VARIABLE 

THKOH 

011400 

REALM 

PROCEDURE 

ROCP 

011362 

REALM 

PROCEDURE 

ASSIGN 

000600 

REALM 

PROCEDURE 

TEMP  10 

311544 

REALM 

VARIABLE 

DENS 

011550 

REALM 

VARIABLE 

Q0 

011554 

REALM 

VARIABLE 

BI. 

01 1560 

REALM 

VARIABLE 

AK 

011564 

REALM 

VARIABLE 

BOIL 

011570 

REALM 

VARIABLE 

ABSORB 

011574 

REALM 

VARIABLE 

JINC 

011600 

INTEGERS 

VARIABLE 

TEMPB 

01  1602 

REALM 

VARIABLE 

ETIME 

011606 

REALM 

VARIABLE 

PCWATR 

011612 

REALM 

VARIABLE 

BLOOD 

011616 

REALM 

VARIABLE 

J 

011622 

INTEGERS 

VARIABLE 
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FORTRAN 

V/ 

STORAGE  MAP 

NAME 

OFFSET 

ATTRIBUTES 

PL2 

011624 

REAL#4 

TRIABLE 

PLN2 

011630 

REAL#4 

TRIABLE 

PL1 

011634 

REAL#4 

TRIABLE 

PLN1 

011640 

REAL#4 

VARIABLE 

DE2 

011644 

REAL*4 

TRIABLE 

DEI 

011650 

REAL#4 

VARIABLE 

I 

011654 

INTEGERS 

TRIABLE 

CLOSE 

000000 

REAL*4 

PROCEDURE 

NFLX 

011656 

INTEGERS 

VARIABLE 

PPL1 

011660 

REAL*4 

VARIABLE 

PPLN1 

011664 

REAL*4 

TRIABLE 

DDE1 

011670 

REAL*4 

VARIABLE 

NXTRA 

011674 

INTEGERS 

TRIABLE 

NXTRA0 

011676 

INTEGERS 

TRIABLE 

ANSR 

011700 

REAL*4 

TRIABLE 

IANSR 

011704 

INTEGERS 

VARIABLE 

TDELT 

011706 

REAL*4 

VARIABLE 

DXTRA 

011712 

REAL#4 

TRIABLE 

ALOG 

000000 

REAL *4 

PROCEDURE 

INUM 

011716 

INTEGERS 

TRIABLE 

AJ 

011720 

REAL*4 

VARIABLE 

Q1 

011724 

SEAL*4 

VARIABLE 

HI 

011730 

REAL#4 

VARIABLE 

DTJ 

011734 

REAL*4 

VARIABLE 

JJJJ 

011740 

INTEGERS 

VARIABLE 

ITTR 

011742 

INTEGERS 

!7AR  TABLE 

IFLX 

011744 

INTEGERS 

VARIABLE 

EITIM1 

011746 

REAL*4 

VARIABLE 

FFDG 

011752 

REAL*4 

VARIABLE 

KFDG 

011756 

INTEGERS 

VARIABLE 

TMPMAX 

011760 

REAL*4 

VARIABLE 

QCONST 

011764 

REAL*4 

VARIABLE 

BLUD 

011770 

REAL*4 

VARIABLE 

M 

011774 

INTEGERS 

VARIABLE 

TIME 

011776 

REAL*4 

TRIABLE 

ITFLG 

012002 

INTEGERS 

I7ARIABLE 

MOD 

000000 

INTEGERS 

PROCEDURE 

P 

012004 

REAL*4 

VARIABLE 

QT 

012010 

REAl*4 

VARIABLE 

N 

012014 

INTEGERS 

VARIABLE 

DT 

012016 

REAL*4 

VARIABLE 

JM1 

012022 

INTEGERS 

VARIABLE 

KK 

012024 

INTEGERS 

TRIABLE 

KHJ 

012026 

INTEGERS 

VARIABLE 

ABS 

000000 

REAL*4 

PROCEDURE 

DW1N 

012030 

REAL*4 

VARIABLE 

EXP 

000000 

REAL*4 

PROCEDURE 

EMTIME 

012034 

REAL*4 

VARIABLE 

AINT 

000000 

REAL*4 

PROCEDURE 

AMOD 

000000 

REAL*4 

PROCEDURE 

EM 

012040 

REAL*4 

VARIABLE 

NN 

012044 

INTEGERS 

VARIABLE 

JLTl 

012046 

INTEGERS 

TRIABLE 

K 

012050 

INTEGERS 

VARIABLE 

DEPTH 

000000 

REAL*4 

PROCEDURE 

TD 

012052 

REAL*4 

TRIABLE 

FORTRAN  Ik  STORAGE  MAP 

NAME  OFFSET  ATTRIBUTES 

I  ERR  012056  INTEGERS  VARIABLE 
IOER  012060  INTEGERS  kARIABLE 
DELEET  000000  REAL*4  PROCEDURE 
EXIT  000000  R£AL*4  ,  PROCEDURE 


FORTRAN  *  11/  1^0 1C-03D  FRI  20-APR-79  10:21:01  •  PAGE  001 

CORE-13K,  UIC-C 120.  73  MODEL3 . MODEL3 -MODEL3 . FTN ; 3 

0001  SUBROUTINE  DEPTH (X.V.N.TD. IERR) 

C 

C  IN17ERSE  INTERPOLATION  ON  TWO  OR  THREE  POINTS  TO  DETERMINE 

C  THRESHOLD  DEPTH  (PREDICTED  BURN  DEPTH)  USING  EITHER 

C  V  OR  LOG(V) 

C 

0002  DIMENSION  X( 1) . V< 1) . 2(3) 

C 

0003  IERR  -  0 

0004  IF  (N.LT.2)  GO  TO  999 

0006  DO  100  1-1. N 

0007  100  Z(I)  -  Y(I) 

0008  Z0  -  1. 

0009  IF  (Z(1).EG.0..OR.Z(2).EQ.0.)  GOTO  140  ! USE  LOGARITHMS? 

0011  IF  (N.£Q.3.AND.Z(3).EQ.0.)  GO  TO  140 

0013  IF  (N.EQ.2)  GO  TO  140 

0015  Z0  -  0.  IUSE  LOGARITHMS 

0016  DO  120  I-l.N 

001?  120  2C I )  -  ALOG(Z(D) 

0018  140  H0  «  Z(2)-Z(l) 

0019  IF  (H0.EQ.0.)  GO  TO  999 

0021  IF  (N.EQ.2)  GO  TO  160 

0023  HI  -  2(3)-Z(2) 

0024  IF  (H1.EQ.0.)  GO  TO  999 

0026  H2  -  Z(3)-Z( 1) 

002?  IF  (H2.EQ.0. )  GO  TO  999 

0029  DZ3  -  Z0-Z(3) 

0030  160  DZ2  -  20-2(2) 

0031  DZ1  •  Z0-Z(1) 

0032  IF  (N.EQ.2)  GO  TO  180 

0034  TD  •  DZl*DZ2*X(3)/(Hi*H2)-DZl*X(2)*DZ3/(N0*Hl)+X( 1)*DZ2*DZ3/ 

1 (H8*H2) 

0035  GO  TO  200 

0036  180  TD  -  (DZ1*X(2)-X(1)*DZ2)/H0 

003?  200  TD  -  EXP(TD) 

0038  RETURN 

0039  999  IERR  -  -l 

0048  RETURN 

0041  END 
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FORTRAN 

W 

STORAGE  MAP 

NAME 

OFFSET 

ATTRIBUTES 

X 

000014 

REAL*4 

PARAMETER 

ARRAV 

(1) 

V 

000016 

SEAL*4 

PARAMETER 

ARRAY 

(1) 

z 

000026 

REAL*4 

ARRAV  (3) 

N 

000020 

INTEGERS 

PARAMETER 

VARIABLE 

TD 

000022 

R£AL#4 

PARAMETER 

TRIABLE 

I  ESS 

000024 

INTEGERS 

PARAMETER 

TRIABLE 

I 

000042 

INTEGERS 

TRIABLE 

Z0 

000044 

REAL#4 

TRIABLE 

ALOG 

000000 

REAL*4 

PROCEDURE 

H0 

000050 

REAL*4 

TRIABLE 

HI 

000054 

REAL*4 

VARIABLE 

H2 

000060 

REAL*4 

TRIABLE 

DZ3 

000064 

REAL*4 

TRIABLE 

DZ2 

000070 

REAL*4 

TRIABLE 

DZ1 

000074 

REAL*4 

TRIABLE 

EXP 

000000 

REAL*4 

PROCEDURE 
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